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Chapter  I:  Introduction 


I  remember  one  beautiful  spring  day  I  was  sitting  in  my  undergraduate  Quantum 
Mechanics  and  Atomic  Physics  class  wishing  that  somehow  we  could  hold  class  outside. 
Of  course  that  would  have  been  impossible  because  this  particular  physics  professor 
always  used  the  blackboard  to  derive  everything  from  first  principles.  We’d  solved 
Schrbdinger’s  equation  for  the  hydrogen  atom  about  a  month  earlier,  and  all  the  lectures 
and  homework  problems  since  had  been  working  toward  solutions  to  multi-electron  and 
multi-nuclear  equations  (where  multi  is  generally  taken  to  mean  two).  Even  taking 
advantage  of  some  minor  approximations  that  the  professor  would  allow  us,  these  were 
long  involved  derivations  filling  reams  of  notebook  paper  and  countless  blackboards  with 
equations.  Perhaps,  on  this  particular  day,  the  weather  was  getting  to  him  as  well,  because 
with  at  least  five  minutes  left  in  the  class  period  he  stepped  back  from  the  board  and  said 
“Well,  you  can  see  where  this  is  going,  and  that  it  is  possible  to  solve  this  case  if  you’re 
willing  to  work  for  it.  But  you  can  also  see  the  situation  gets  complex  rather  quickly. 
What  really  impresses  me  is  that  the  guys  across  the  street  (in  the  chemistry  building)  can 
solve  this  for  several  atoms  and  dozens  of  electrons.  Don’t  ask  me  how  though.  It’s  left 
as  an  exercise  for  the  inspired  student  to  go  ask  them.” 

Well,  I  didn’t  walk  across  the  street  and  ask.  Instead  I  called  up  my  friend.  Brad 
DeLeeuw,  who  was  then  a  first  year  graduate  student  in  Quantum  Chemistry.  He  did  his 
best  to  explain  everything  to  me  over  the  phone,  but  what  I  got  out  of  it  was  that  Quantum 
Chemists  took  advantage  of  a  couple  more  useful  approximations  and  then  used  computers 
to  do  the  bulk  of  the  work.  The  thought  of  letting  computers  do  the  work  for  me  had  a  nice 
ring  to  it,  so  to  learn  more  I  joined  the  summer  undergraduate  research  program  at  the 
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Center  for  Computational  Quantum  Chemistry.  In  my  naivete  I  decided  to  have  the 
computers  solve  the  entire  C2H2O  potential  energy  hypersurface  and  determine,  among 
other  things,  if  oxirene  was  the  symmetric  stable  intermediate  in  the  Wolff  decomposition 
of  diazoketone  to  ketene.  Over  the  course  of  the  summer  I  barely  had  time  to  uncover  the 
true  difficulties  of  that  surface  and  didn’t  have  a  real  chance  to  solve  any  of  them.  And  the 
problems  of  that  surface  have  continually  haunted  me  for  most  of  my  graduate  career. 

Dr.  Yamaguchi  tried  to  help  me  to  reach  some  solutions  in  the  oxirene  problem 
during  my  first  year  of  graduate  school.  He  too  decided  that  the  problems  of  oxirene  were 
too  difficult,  so  he  turned  my  attention  to  an  easier  set  of  species  for  which  I  might  expect 
to  successfully  complete  a  project,  the  triplet  acetylene  and  vinylidene  surfaces.  Along  with 
ammonia  and  formaldehyde,  acetylene  is  one  of  the  three  tetraatomic  molecules  most 
exhaustively  characterized  by  spectroscopic  techniques.  Interestingly,  though,  many 
important  energetic  and  spectral  concerns  for  the  low  lying  excited  states  of  acetylene  and 
their  vinylidene  isomers  were  not  yet  determined  in  the  laboratory.  So  we  undertook  a 
study  to  provide  as  much  theoretical  insight  as  possible  into  that  area  of  chemistry.  The 
results  of  that  research  are  presented  in  Chapter  H. 

During  my  second  year  of  graduate  school,  another  graduate  student,  Jeff  Pilgrim, 
approached  me  with  a  question.  In  a  study  of  aluminum  and  water  clusters,  he  and 
Professor  Duncan  had  uncovered  an  unexpected  spectrum  at  44  amu.  Undoubtedly  it  was 
either  HAIO  or  AlOH,  but  which?  There  wasn’t  enough  available  data  in  any  journal 
articles  for  them  to  form  a  decision.  Jeff  asked  if  it  were  possible  for  me  to  quick  calculate 
what  he  needed  to  know.  I  assured  him  that  I  would  and  sent  him  on  his  way.  I  got  him 
those  results  (though  perhaps  not  as  quickly  as  he  envisioned)  and  they  are  presented  here 
in  Chapter  III. 

With  those  research  projects  successfully  under  my  belt,  I  returned  my  focus  to  the 
evils  of  the  question  about  oxirene  participation  in  the  ketene  rearrangement.  Although  I 
did  wish  to  see  the  logic  of  science  triumph  over  the  eomplexities  of  nature,  my  real 
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motives  for  returning  to  this  problem  were  not  completely  altruistic.  Rather,  everyone  who 
knew  I  had  started  the  oxirene  project  (not  the  least  of  whom  was  Professor  Schaefer) 
wanted  to  know  the  results.  The  worst  aspect  of  the  study  of  the  entire  ketene  -  oxirene 
isomerization  is  that  different  levels  of  quantum  chemical  theory  give  different  results. 
Some  structures  that  are  stationary  points  on  the  potential  energy  hypersurface  at  one  level 
of  theory  disappear  at  other  levels  of  theory.  The  oxirene  stationary  point  itself  switches 
from  an  equilibrium  to  a  transition  state  with  changes  in  the  level  of  theory.  Determination 
of  which  level  of  theory  gives  the  correct  qualitative  results  was  a  major  problem. 
Presented  in  Chapter  IV  is  my  study  of  the  oxirene  potential  energy  surface  to  decide  the 
best  level  of  theory  at  which  to  study  the  entire  isomerization  reaction  surface. 


Chapter  II 


Isomerization  Reactions  on  the  Lowest 
Potential  Energy  Hypersurface  of  Triplet 
Vinylidene  and  Triplet  Acetylene^ 


1  G.  V.  Vacek,  J.  R.  Thomas,  B.  J.  DeLeeuw,  Y.  Yamaguchi  and  H.  F.  Schaefer  7.  Chem.  Phys.  98, 
4766  (1993). 
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Triplet  vinylidene,  first  predicted  to  have  a  sizable  barrier  to 
unimolecular  rearrangement  in  1978  by  theory,  has  now  been  observed 
under  three  different  sets  of  experimental  conditions.  In  order  to 
quantitatively  characterize  the  potential  energy  hypersurface  of  triplet 
vinylidene  and  triplet  acetylene,  high-level  ab  initio  quantum  mechanical 
methods  have  been  employed.  Basis  sets  as  large  as  triple  zeta  plus  two 
sets  of  polarization  functions  augmented  with  higher  angular  momentum 
functions  [TZ(2df,2pd)]  have  been  utilized  in  conjunction  with  correlated 
methods  as  sophisticated  as  the  coupled  cluster  approach  including  all 
single,  double  and  perturbative  triple  excitations  [CCSD(T)].  Of  particular 
interest  are  predictions  of  the  zero-point  vibrational  energy  corrected  barriers 
for  rearrangement  of  a  ^82  vinylidene  to  b  ^Bu  trans-bent  acetylene  and  of 
cis-bent  a  ^B2  acetylene  to  trans-bent  b  ^Bu  acetylene.  At  the  highest  level 
of  theory  used  here,  TZ(2df,2pd)  CCSD(T),  these  are  predicted  to  be  47.9 
kcal/mole  and  13.0  kcal/mole,  respectively.  The  physical  properties  of 
a  3B2  vinylidene,  including  the  dipole  moment,  harmonic  vibrational 
frequencies  and  the  infrared  intensities,  have  also  been  reported. 


Introduction 

For  nearly  a  decade  it  was  unclear  whether  singlet  vinylidene 


C  : 


(1) 


was  a  shallow  minimuml’^  on  the  C2H2  ground  state  potential  energy  hypersurface  or 
merely  a  transition  state4>5  for  the  degenerate  rearrangement  of  hydrogen  atoms,  e.g. 


H_c12_q13 - D  ^ 


H 

\ 

C'2=Z=C'3 


D— — H 


(2). 
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Recent  negative  ion  photoelectron  spectroscopy  experiments  by  Ervin,  Ho  and  Lineberger^ 
have  shown  that  the  vinylidene  potential  well  must  be  at  least  as  deep  as  the  experimentally 
observed  2<—0  CH2  rock  transition,  namely  450  cm'*  or  1.3  kcal/mole.  The  most  recent 
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and  extensive  theoretical  study  by  Gallo,  Hamilton,  and  Schaefer^  predicted  a  barrier  of  1.6 
kcal/mole  including  corrections  for  zero-point  vibrational  energies  (ZPVE's). 

In  the  case  of  triplet  states  of  vinylidene,  the  situation  is  still  somewhat  nebulous. 
The  first  solid  evidence  that  the  a  ^62  state  of  vinylidene  might  be  experimentally 
observable  was  provided  by  theory.  In  1978,  Conrad  and  Schaefer^  predicted  a  classical 
barrier  of  56  kcal/mole  for  its  isomerization  to  the  appropriate  triplet  state  of  acetylene  using 
a  basis  set  and  method  for  the  treatment  of  electron  correlation  effects  that  are  of  moderate 
quality  by  current  standards.  They  also  estimated  the  activation  energy  (including  ZPVE 
corrections)  to  be  about  50  kcal/mole.  Shortly  thereafter  in  flash  photolysis  experiments, 
Laufer^  identified  a  reactive  intermediate  with  composition  C2H2  as  the  a  ^62  state  of 
vinylidene.  Laufer  and  coworkers^®  subsequently  characterized  the  reactivity  and 
quenching  rates  of  this  triplet  state  of  vinylidene.  In  related  work,  Siilzle  and  Schwarz^^ 
observed  a  triplet  state  of  vinylidene  with  a  lifetime  greater  than  0.4  (xsec  by  neutralization- 
reionization  mass  spectrometry. 

The  most  quantitative  energetic  information  concerning  triplet  vinylidene  was 
provided  by  the  negative  ion  photoelectron  spectroscopy  studies  of  Lineberger  and 
coworkers. 6.12  They  determined^  singlet-triplet  splittings  of  To(a  ^62)  =  47.6  ±  0.1 
kcal/mole  and  To(b  ^A2)  =  63.5  ±  0.5  kcal/mole  for  the  two  lowest  excited  states  of 
vinylidene.  Thus  their  experimentally  observed  value  of  ATo(b  3A2  -  a  3B2)  =  15.9  ±  0.6 
kcal/mole  was  obtained,  in  reasonable  agreement  with  the  earlier  theoretical  predictions  of 
ATe(b  3A2  -  a  3B2)  =  15.2  kcal/mole  by  Dykstra  and  Schaeferi  and  ATo(b  ^A2  -  a  3B2)  = 
13.9  kcal/mole  by  Osamura  and  Schaefer.  13  The  very  earliest  prediction  of  the  3B2  -  3A2 
energy  separation  by  Davis,  Goddard  and  Hardingl'l  was  26.6  kcal/mole. 

In  1981  an  interesting  qualitative  molecular  orbital  discussion  of  triplet  vinylidene 
was  presented  by  Harding  13  in  the  context  of  his  excellent  study  of  the  vinyl  and  ethyl 
radicals  and  triplet  methylcarbene.  Noting  that  the  doubly-occupied  5ai  lone  pair  orbital 
and  an  unoccupied  p-7t  orbital  in  the  b  3A2  state  of  vinylidene  have  characteristics  which 
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strongly  resemble  the  doubly-occupied  3ai  orbital  and  nonbonding  p-type  orbital  in  singlet 
methylene,  he  concluded  that  the  b  3A2  state  "should  possess  a  low  barrier  to  hydrogen 
migration  analogous  to  that  of  other  singlet  carbenes."  Harding  also  claimed  that  the 
transition  state  for  isomerization  to  acetylene  on  the  3A2  surface  should  be  highly 
nonplanar.  This  leads  to  an  avoided  crossing  of  the  ^62  and  potential  energy  surfaces 
(PES’s)  at  the  Ci-symmetry  transition  state  for  1,2-hydrogen  shift. 

Interaction  between  theory  and  experiment  provided  abundant  evidence  for  the 
existence  of  low-lying  triplet  excited  electronic  states  of  acetylene.  For  example,  theoretical 
predictions^  followed  quickly  by  experimental  confirmationS7  left  Hftie  doubt  that  a  cw-bent 
^B2  state  is  the  lowest  excited  state  of  acetylene.  Nevertheless,  Dupre,  lost,  Lombardi, 
Green,  Abramson  and  Fields^  have  written  the  only  experimental  paper  addressing  these 
triplet  states  during  the  past  three  years.  Surprisingly,  the  singlet-triplet  splitting  To(a  3B2) 
has  not  yet  been  determined  experimentally  for  acetylene.  In  addition,  no  experimentally 
observed  spectroscopic  features  have  been  assigned  to  other  low-lying  electronic  states, 
such  as  the  predicted^^  trans-bent  b  ^Bu  state.  Theoretical  work  by  Lischka  and  Karpfen^^ 
in  1986  and,  more  recently,  by  Yamaguchi,  Vacek,  and  Schaefer^O  constitute  the  only 
high-level  studies  concentrating  on  these  states  since  the  ground-breaking  1978  theoretical 
study.  Both  of  these  later  works  have  focused  on  the  molecular  properties  of  stable 
minima  on  the  triplet  acetylene  PES. 

Based  upon  a  recent  analysis  of  Zeeman  anticrossing  (ZAC)  spectra  of  gas  phase 
acetylene,  it  has  been  hypothesized  that  the  Si  excited  singlet  state  of  acetylene  must  lie 
close  in  energy  to  at  least  one  of  the  isomerization  transition  states  on  the  triplet  PES  of 
acetylene.  Energetically  speaking,  the  two  most  likely  candidates  are  the  transition  state  for 
cw-bent  a  ^B2  to  trans-bent  b  ^Bu  acetylene  isomerization  and  the  transition  state  for 
1,2-hydrogen  migration  between  a  3B2  vinylidene  and  b  ^Bu  acetylene.  Information  gained 
from  a  high-level  theoretical  study  of  these  two  transition  states  and  the  a  ^B2  state  of 
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vinylidene  should  prove  quite  useful  in  further  experimental  work  on  excited  states  of 
acetylene. 

In  light  of  the  large  body  of  experimental  results, ^>9  it  is  remarkable  that  the  1978 
work  by  Conrad  and  Schaefer^  is  the  only  theoretical  study  of  1,2-hydrogen  migration  in 
triplet  vinylidene.  Although  theory  assisted  Laufer's  1980  discovery^  of  a  ^62  vinylidene, 
it  is  apparent  that  experiment  has  now  taken  a  clear  lead.  Here  we  attempt  to  bring  theory 
back  as  a  synergetic  partner.  Ervin  et  al.^  have  observed  three  (out  of  six  possible) 
fundamental  vibrational  frequencies  for  a  ^62  vinylidene.  All  six  harmonic  vibrational 
frequencies  are  reliably  predicted  here.  Finally,  quantitative  predictions  of  the  activation 
energy  for  1,2-hydrogen  migration  in  a  ^62  vinylidene  and  the  related  barrier  for  cis-trans 
triplet  acetylene  interconversion  are  made. 

Electronic  Structure  Considerations 

Ground  state  vinylidene  is  described  in  C2v  symmetry  by  the  electronic 
configuration 

[core](3ai)2(4ai)2(lb2)2(5ai)2(lbi)2  X  ^Ai  vinylidene  (3). 

Note  that  [core]  abbreviates  the  two  lower-lying  C  Is-like  orbitals.  The  CC  double  bond  in 
this  structure  consists  of  the  3ai  a-bonding  (a*^)  and  the  Ibi  7C-bonding  (tt^)  orbitals.  The 
highest  occupied  molecular  orbital  (HOMO)  is  the  Ibi  7i*^  orbital.  The  5ai  lone  pair  on  the 
terminal  (carbene)  carbon  is  the  second-highest  occupied  molecular  orbital  (SHOMO). 
These  two  HOMOs  are  energetically  very  close.  The  lowest  unoccupied  molecular  orbital 
(LUMO),  the  2b2  orbital,  is  essentially  an  in-plane  nonbonding  p-type  orbital  centered 
upon  the  terminal  carbon.  In  this  particular  case,  the  electronic  configuration  of  the  lowest 
excited  state  of  vinylidene,  the  a  3B2  state,  differs  from  that  of  ground  state  vinylidene  by 
excitation  of  a  single  electron  from  the  SHOMO  to  the  LUMO.  Interestingly,  exciting  an 
electron  from  the  HOMO  to  the  LUMO  of  X  ^Ai  vinylidene  leads  to  the  second-lowest 
b  3A2  excited  state  of  vinylidene.  The  reasons  for  this  unusual  ordering  have  been 
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discussed  adequately  by  Dykstra  and  Schaefer.  1  Written  in  terms  of  C2v-symmetry  orbitals 
the  electronic  configurations  of  these  two  states  are  thus: 

[core](3ai)2(4ai)2(lb2)^(lbi)2(5ai)(2b2)  a  ^62  (4); 

[core](3ai)2(4ai)2(lb2)^(5ai)^(lbi)(2b2)  b3A2  (5). 

Based  upon  the  electronic  configuration  (4)  given  above,  the  a  ^32  state  is  expected  to  have 
a  genuine  CC  double  bond  since  both  the  3ai  and  Ibi  orbitals  remain  doubly-occupied. 
Because  of  the  singly-occupied  it'^  orbital  in  its  electron  configuration  (5),  the  b  3A2  state  is 
best  considered  to  have  a  bond  order  of  one  and  a  half 

In  its  linear  ground  state,  acetylene  can  be  described  by  the  electronic  configuration 
[core](2ag)2(2au)2(3ag)2(ljCu)4  X  iSg'*’  acetylene  (6). 

The  liUu  HOMO  of  this  ground  state  is  a  orbital,  and  the  lllg  LUMO  is  an  antibonding 
(7t*)  orbital.  When  acetylene  is  distorted  to  C2v  symmetry,  the  degenerate  liiu  HOMO 
splits  into  a  lower-energy  Ibi  orbital  and  a  higher-energy  4ai  orbital,  while  the  liCg  LUMO 
becomes  the  3b2  and  la2  orbitals  (with  the  3b2  lower  in  energy).  Thus,  written  in  terms  of 
C2v-symmetry  orbitals,  the  electronic  configuration  of  ground  state  acetylene  is 

[core](2ai)2(2b2)2(3ai)2(lbi)2(4ai)2  (7). 

The  lowest  excited  triplet  state  has  a  cw-bent  structure  and  belongs  to  the  C2v  point  group. 
Its  electronic  configuration  is  derived  from  (7)  by  excitation  of  a  single  electron  from  the 
HOMO  to  the  LUMO: 

[core](2ai)2(2b2)2(3ai)2(lbi)2(4ai)(3b2)  5^32  (8). 

Distortion  of  linear  ground  state  acetylene  to  trans-hcnt  C2h  symmetry  will  cause  a  different 
splitting  of  the  degenerate  orbitals,  namely  7t^  lau,3bu  and  n*  — >  4ag,lbg.  Written  in 
terms  of  C2h-synimetry  MO's,  the  electronic  configuration  becomes 

[core](2ag)2(2bu)2(3ag)2(lau)2(3bu)2  (9) 

for  the  ground  state  of  acetylene.  Excitation  of  an  electron  from  the  3bu  HOMO  to  the  4ag 
LUMO  results  in  the  second-lowest  triplet  state  of  acetylene.  This  triplet  state  has  a 
trans-hQni  stmcture  with  the  electronic  configuration 
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[core](2ag)2(2bu)2(3ag)2(lau)2(3bu)(4ag)  b3Bu  (10). 

The  X  IZg'^  state  of  acetylene  has  a  prototypical  triple  bond  (one  a  bond  and  two  it  bonds). 
In  both  excited  states  of  acetylene  discussed  above,  the  carbon-carbon  bond  involves  two 
electrons  in  the  orbital,  three  in  tt*’  orbitals,  and  one  in  a  7t*  orbital.  Therefore,  the  a  ^62 
and  b  ^Bu  states  of  acetylene  are  expected  to  exhibit  double-bond  characteristics. 

Under  the  constraint  of  Cs  symmetry  with  the  %  symmetry  element  defined  by  the 
molecular  plane,  the  qualitative  electronic  configurations  of  a  3B2  vinylidene,  a  3B2 
acetylene  and  b  ^Bu  acetylene  may  all  be  written  as 

[core](3a')2(4a')2(5a')2(la")2(6a')(7a’)  ^a’  (11). 

Though  they  may  have  identical  qualitative  electronic  configurations,  each  MO  may  have  a 
significantly  different  nature  in  each  of  these  three  triplet  states.  Incidentally,  the  qualitative 
electronic  configuration  of  any  planar  transition  state  between  these  three  triplet  states  on 
the  C2H2  triplet  surface  is  also  represented  by  (11).  On  the  other  hand,  written  in  terms  of 
Cs-symmetry  MO's  with  the  molecular  plane  again  defining  Oh.  the  b  3A2  state  of 
vinylidene  has  the  electronic  configuration 

[core](3a')2(4a’)2(5a')2(6a')2(la’')(7a’)  3a"  (12). 

Distortion  of  the  molecular  electronic  states  described  by  equations  (11)  and  (12) 
into  Cl  symmetry  changes  the  qualitative  electronic  configurations  to 

[core](3a)2(4a)2(5a)2(6a)2(7a)(8a)  3a  (13). 

Thus  all  the  triplet  states  of  C2H2  mentioned  above  become  indistinct  when  distorted  into 
Cl  symmetry. 

Theoretical  Procedures 

Three  basis  sets  were  employed  in  this  study.  The  first  basis  set  was  of  double-zeta 
plus  polarization  (DZP)  quality.  It  consisted  of  a  standard  Huzinaga-Dunning  double-zeta 
get2l,22  of  contracted  Gaussian  functions  [designated  (9s5p/4s2p)  for  carbon  and  (4s/2s) 
for  hydrogen]  augmented  with  a  single  set  of  d-like  polarization  functions  on  C  and  p 


11 


functions  on  H.  The  exponents  of  the  polarization  functions  were  ad(C)  =  0.75  and  ap(H) 
=  0.75.  The  triple-zeta  plus  two  sets  of  polarization  functions  (TZ2P)  basis  set  was 
composed  of  Dunning's^^  (5s3p)  contraction  of  Huzinaga's^l  (10s6p)  primitive  set  of 
Gaussian  functions  for  carbon  and  Dunning's^^  (3s)  contraction  of  Huzinaga's^^  (5s)  set 
for  hydrogen  augmented  with  two  sets  of  polarization  functions  for  each  atom.  The 
polarization  function  orbital  exponents  were  0C(j(C)  =  1.50,  0.375  and  Ctp(H)  =  1.50, 
0.375.  The  final  basis  set  was  formed  by  supplementing  the  TZ2P  basis  with  higher 
angular  momentum  polarization  functions  [af(C)  =  0.8  and  ad(H)  =  1.0].  This  basis  set 
has  been  designated  TZ(2df,2pd).  Throughout  this  study,  sets  of  six  Cartesian  d-like  and 
ten  Cartesian  f-like  polarization  functions  were  employed. 

The  molecular  structures  of  the  lowest  triplet  state  of  vinylidene  and  the  transition 
states  for  isomerization  reactions  were  fully  optimized  using  analytic  gradient  techniques 
for  restricted  Hartree-Fock  self-consistent  field  (SCF)24,25  and  configuration  interaction 
(Cl)  wavefunctions.26.27  in  all  cases,  the  residual  Cartesian  and  internal  coordinate 
gradients  were  less  than  10-6  atomic  units.  Analytic  second  derivative  procedures28.29 
were  used  to  evaluate  the  quadratic  force  constants  for  the  SCF  wavefunctions. 

Cl  and  coupled  cluster  (CC)  methods  were  employed  to  include  the  effects  of 
electron  correlation.  All  single  and  double  excitations  from  the  SCF  reference 
configuration  were  included  (CISD  and  CCSD),  with  the  exception  that  only  the  valence 
electrons  were  explicitly  correlated.  In  this  manner,  the  two  lowest-lying  MOs  (C  Is-like) 
were  held  doubly-occupied  (frozen  cores),  and  the  two  highest- lying  orbitals  (C  ls*-like) 
were  excluded  (deleted  virtuals)  from  the  configuration  space.  With  the  TZ(2df,2pd)  basis 
set,  this  correlation  scheme  involved  38,585  configurations  for  the  C2v-symmetry 
wavefunction  of  ^62  vinylidene,  144,777  configurations  for  the  Ci-symmetry  transition 
state  from  ^B2  vinylidene  to  acetylene,  and  76,557  configurations  for  the  Cs-symmetry 
transition  state  to  interconversion  between  ^B2  acetylene  and  ^Bu  acetylene  in  the  Hartree- 
Fock  interacting  space^O.Bl  for  the  CISD  and  CCSD  wavefunctions.  The  shape  driven 
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graphical  unitary  group  approach^^  ^as  used  to  obtain  CISD  wavefunctions.  An  estimate 
of  the  contribution  of  unlinked-cluster-quadruple-excitations  to  the  CISD  energies  was 
obtained  via  Davidson's  correction.^^’^'^  Such  results  are  designated  CISD+Q.  The 
method  of  finite  central  differences  of  analytic  gradients26.27  ^^as  utilized  to  evaluate  the 
quadratic  force  constants  for  the  CISD  wavefunctions. 

Open-shell  CCSD  wavefunctions  were  constructed  using  the  methodology 
developed  by  Scuseria.35  jhis  follows  an  approach  originally  developed  by  Purvis  and 
Bartlett36  for  the  closed-shell  case  and  modified  by  Scuseria  and  Janssen.37  xhe  effect  of 
connected  triple  excitations  was  included  using  the  perturbative  approach  suggested  by 
Raghavachari,  Trucks,  Pople  and  Head-Gordon  [CCSD(T)]38  and  implemented  in  our 
group  by  Scuseria.35  Improved  predictions  of  relative  energies  with  a  particular  basis  set 
were  obtained  by  evaluating  coupled  cluster  single-point  energies  at  the  corresponding 
CISD  optimized  geometries.  The  textual  notation  used  to  designate  such  single-point 
energy  calculations  in  this  paper  is  as  follows:  a  DZP  CCSD(T)  single-point  energy 
evaluated  at  the  DZP  CISD  optimized  geometry  will  be  referred  to  as  DZP 
CCSD(T)//CISD. 

Results 

Equilibrium  structures  predicted  for  a  ^62  vinylidene  are  presented  in  Figure  1. 
The  optimized  geometries  for  the  two  "true"  transition  state  structures  (i.e.  those  having 
only  one  imaginary  frequency)  are  shown  in  Figures  2,  3  and  4.  Two  other  stationary 
points  that  were  examined  on  the  DZP  SCF  surface  and  found  to  have  more  than  one 
imaginary  frequency  are  discussed  below  but  not  depicted. 

The  harmonic  vibrational  frequencies  predicted  for  the  three  stationary  points,  as 
well  as  the  total  energies,  dipole  moments  and  ZPVEs  are  listed  in  Tables  I,  III  and  IV. 
The  corresponding  infrared  (IR)  intensities  for  a  ^62  vinylidene  are  also  given  in  Table  I. 
Harmonic  vibrational  frequencies  for  isotopically  substituted  a  ^62  di-  and  d2-vinylidenes 
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are  presented  in  Table  H.  It  should  be  noted  that  the  vibrational  frequencies  predicted  by 
this  research  are  harmonic  frequencies  in  all  cases,  not  fundamentals;  thus  the  contribution 
of  anharmonic  terms  is  not  included  in  our  results. 

The  CISD+Q,  CCSD,  and  CCSD(T)  total  energies  determined  at  CISD  optimized 
geometries  are  shown  in  Table  V.  Relative  energies  of  a  3B2  vinylidene  with  respect  to 
X  ^Ai  vinylidene  given  in  Table  VI  were  obtained  using  energetic  information  from  Tables 
I  and  V  as  well  as  results  published  by  Gallo  et  alP  where  possible.  Relative  energies  of 
the  other  studied  stationary  points  with  respect  to  a  ^62  acetylene  are  summarized  in  Table 
Vn  and  have  been  obtained  using  the  energetic  information  presented  in  Tables  I,  and 
in  -  V.  Relative  energies  of  the  b  ^Bu  state  of  acetylene  were  also  included  in  Table  VH  to 
allow  more  detailed  comparison.  The  values  used  for  these  two  triplet  states  of  acetylene 
are  taken  directly  from  our  previous  study^O  for  identical  basis  sets  and  methods.  For  the 
CISD+Q,  CCSD  and  CCSD(T)  energies,  the  CISD  ZPVE  values  for  the  corresponding 
basis  sets  were  employed  to  make  best-estimate  corrections. 

In  some  cases,  no  previously  reported  information  was  available  for  X  ^Ai 
vinylidene,  a  ^B2  acetylene  and  b  ^Bu  acetylene  at  levels  of  theory  consistent  with  this 
work.  For  these  states,  previously  unreported  harmonic  vibrational  frequencies  and 
energetic  data  determined  in  this  research  are  presented  in  Appendix  A.  None  of  these 
results  were  qualitatively  different  from  those  reported  in  earlier  works^’20  for  the  relevant 
states. 

Discussion 

Geometries 

d  ^82  Vinylidene 

Based  upon  the  results  of  their  Franck-Condon  simulations  of  observed  spectra, 
Ervin  et  al.^  reported  an  "experimentally  determined"  molecular  geometry  for  a  '^82 
vinylidene:  re(CC)  =  1.346  +  0.040  A,  re(CH)  =  1.090  ±  0.009  A,  and  OelHCH)  =  1 18.9 
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±  2.7°.  The  Franck-Condon  simulations  were  optimized  relative  to  the  molecular  structure 
of  X  lAi  vinylidene  predicted^  at  the  TZ  +  P  CISD  level  of  theory.  Therefore,  it  is  not 
very  surprising  that  all  of  the  CISD  geometric  parameters  for  the  a  ^62  state  of  vinylidene 
shown  in  Figure  1  fall  safely  within  the  experimentally  reported  error  margins. 

A  comparison  of  the  a  3B2  vinylidene  geometry  in  Figure  1  with  that  of  X  ^Ai 
vinylidene  reported  by  Gallo  et  alP  at  comparable  levels  of  theory  reveals  the  small  changes 
in  geometry  that  result  from  exciting  an  electron  from  the  5ai  jp-type  lone  pair  MO  to  the 
2b2  nonbonding  p-type  MO.  The  CC  bond  slightly  elongates  in  a  3B2  vinylidene  [re(CC)  = 
1.3121  A  at  TZ2P  CISD]  relative  to  the  ground  state  [re(CC)  =  1.2957  A  at  TZ2P  CISD],^ 
but  both  states  obviously  involve  CC  double  bonds.  (In  ethylene,  the  experimental  double¬ 
bond  length  is  ro(CC)  =  1.337  The  HCH  angle  is  2°  smaller  in  the  a  3B2  state  than  in 
the  X  lAi  state  due  to  interaction  between  the  singly-occupied  2b2  orbital  and  the  lb2  CH 
orbital.  There  is  no  significant  change  in  the  CH  bond  lengths  between  the  a  ^^2  and  X 
lAi  states  of  vinylidene. 

Transition  State  for  Triplet  Vinylidene  Triplet  Acetylene  Isomerization 

An  initial  search  was  made  for  a  planar  transition  state  for  the  ^B2  vinylidene  4-> 
3Bu  acetylene  isomerization  reaction.  A  stationary  point  with  a  "bridged"  structure  similar 
to  that  previously  reported  by  Conrad  et  alf  was  readily  found.  At  the  DZP  SCF  level,  the 
optimized  geometry  is  rc(CiC2)  =  1.291  A,  re(CiHbridged)  =  1.250A,  re(C2Hnonbridged)  = 
1.089A,  ee(C2CiHbridged)  =  60.5°,  0e(CiC2Hnonbridged)  =  135.3°.  Vibrational  analysis  of 
this  structure,  however,  revealed  a  2780i  cm'l  a'  imaginary  mode  corresponding  to 
in-plane  motion  of  the  migrating  (bridging)  hydrogen  and  a  735i  cm'l  a"  imaginary  mode 
corresponding  to  an  HCCH  torsional  motion.  This  Cj-symmetry  structure  is  thus  properly 
characterized  as  a  higher-order  saddle  point  with  a  Hessian  index  of  two.  The  true 
isomerization  transition  state  with  a  Hessian  index  of  one  should  therefore  belong  to  the 
Ci-symmetry  point  group. 
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When  the  Cs-symmetry  constraint  was  removed,  the  genuine  Ci -symmetry 
isomerization  transition  state  was  located.  Its  structure  is  shown  in  Figure  2.  The  CC 
bond  length  predicted  for  this  structure  (e.g.  1.410A  for  TZ2P  CISD)  implies  that  the 
transition  state  structure  involves  a  bond  of  order  one  and  a  half  between  the  carbon  atoms. 
This  is  distinctly  different  from  the  doubly-bonded  [re(CC)  =  1.244A  at  TZ2P  CISD]”^ 
Cs-symmetry  transition  state  structure  predicted  for  isomerization  of  X  vinylidene  to 
X  acetylene.  Another  notable  difference  between  the  C2H2  singlet  ground-state  PES 
and  this  excited  triplet  state  PES  is  that  the  transition  state  structure  for  1,2-hydrogen 
migration  in  the  former  is  planar^  while  it  is  nonplanar  in  the  latter.  In  fact,  the  torsional 
angle  between  the  HCC  planes  for  the  latter  (depicted  in  Figure  3)  is  115.7°  at  the 
TZ(2df,2pd)  CISD  level  of  theory.  Thus,  in  the  case  of  the  1,2-hydrogen  shift  transition 
state  for  a  ^62  vinyhdene,  the  migrating  hydrogen  lies  in  a  position  nearly  perpendicular  to 
the  plane  of  the  other  three  atoms  rather  than  lying  in  that  plane. 

The  geometrical  differences  between  the  isomerization  transition  state  structures  on 
the  singlet  ground-state  PES  and  the  triplet  excited-state  PES  can  be  attributed  to  two 
effects.  First,  in  the  nonplanar  triplet  isomerization  transition  state,  the  CC  Jt  bond  breaks 
and  the  methylene  carbon  hybridization  changes  as  the  CHC  two-electron-three-center  bond 
forms  during  1,2  hydrogen  migration.  Second,  as  Harding  stated  so  succinctly,  there  is 
an  avoided  crossing  of  the  two  lowest  triplet  state  PES's  in  Ci -symmetry.  The 
1,2-hydrogen  shift  transition  state  for  b  3A2  vinylidene  should  have  a  long  CC  single  bond 
and  an  overall  structure  that  is  strongly  nonplanar.  The  mixing  of  these  two  states 
produces  an  isomerization  transition  state  for  the  lowest  triplet  state  with  a  longer  CC  bond 
distance  and  a  nonplanar  structure.  For  the  lowest  triplet  state  surface,  this  mixing  also 
significantly  lowers  the  activation  barrier  to  hydrogen  migration  via  a  Ci-symmetry  path 
relative  to  that  of  the  Cs-symmetry  reaction  pathway. 


16 


Transition  State  for  cis-trans  Triplet  Acetylene  Isomerization 

Initial  searches  for  the  isomerization  transition  state  between  b  ^Bu  and  a  ^62 
acetylene  were  performed  along  two  paths  of  inquiry.  The  first  cis-trans  isomerization 
transition  state  was  sought  under  C2-symmetry  constraints  with  interconversion  taking 
place  through  rotation  about  the  CC  axis.  Another  possible  isomerization  path  was  an 
in-plane  inversion.  The  second  search  was  therefore  made  under  the  constraint  of 
Cs-symmetry. 

The  search  for  a  C2-symmetry  transitions  state  at  the  DZP  SCF  level  led  to  a 
structure  with  linear  geometry:  re(CC)  =  1.327A  and  re(CH)  =  l.OSSA.  Vibrational 
analysis  revealed  that  this  structure  has  a  Hessian  index  of  two,  indicating  that  it  is  a 
higher-order  saddlepoint  and  not  a  true  transition  state.  The  first  imaginary  frequency,  a 
trani'-bending  a-symmetry  mode,  had  a  magnitude  of  1379i  cm'l.  The  second  imaginary 
frequency,  1550i  cm"l,  corresponds  to  a  cw-bending  b-symmetry  mode. 

The  true  triplet  cis-trans  isomerization  transition  state,  illustrated  in  Figure  4,  was 
found  to  have  Cg  symmetry.  The  geometry  of  this  structure  is  about  what  one  would  expect 
for  an  inversion  isomerization  reaction.  The  inverting  hydrogen  is  very  nearly  collinear 
with  the  CC  bond  axis,  but  it  is  4-6°  rrans-bent  with  respect  to  the  other  hydrogen.  This 
structure  agrees  with  Hammond's  Postulate'*®  since  the  cis-hent  ^B2  —>  trans-bent  ^Bu 
acetylene  interconversion  is  an  endothermic  reaction.  The  bond  lengths  and  angles  are 
reasonably  close  to  the  related  geometrical  parameters  of  both  the  cis-  and  trans-bent 
minima,20  with  the  exception  of  the  CH  bond  length  of  the  inverting  hydrogen.  In  the 
a  3B2  and  b  ^Bu  triplet  states  of  acetylene,  the  orbitals  of  the  carbon  atoms  in  the  CH  bonds 
may  be  thought  of  as  typical  j'p^.hybridized  orbitals.  In  the  inversion  transition  state 
structure,  however,  the  orbital  of  the  carbon  atom  in  the  CHinverting  bond  is  best  viewed  as 
an  sp-hybrid  orbital.  The  greater  5-character  leads  to  the  shorter  CHinverting  bond  length 
relative  to  the  other  CH  bond  length. 
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Harmonic  Vibrational  Frequencies 

d  ^B2  Vinylidene 

The  three  experimentally  observed  ai  fundamental  frequencies  reported  by  Ervin 
et  al.^  [CH  symmetric  stretch  2930  ±  10  cm'^,  CC  stretch  1530  ±  70  cm'^  and  CH2 
scissors  1375  ±  10  cm'l]  agree  well  with  the  corresponding  harmonic  frequencies  listed  in 
Table  I.  For  example,  at  the  TZ(2df,2pd)  CISD  level  of  theory  we  have  predicted  the 
harmonic  frequencies  for  the  corresponding  ai  modes  to  be  3143  cm’^  1637  cm-l  and 
1451  cm'l.  If  a  harmonic-to-fundamental  scaling  factor  is  employed  (typically  0.92  -  0.96 
at  levels  including  electron  correlation), the  agreement  becomes  remarkably  good.  In 
the  case  of  a  3B2  C2D2  (see  Table  II),  Ervin  et  al.^  reported  fundamental  frequencies  of 
2160  ±  10  cm-1, 1495  ±  10  cm-^  and  1010  ±  10  cm'l  for  the  CD  symmetric  stretch,  the  CC 
stretch  and  CD2  scissors  modes,  respectively.  At  the  TZ(2df,2pd)  CISD  level  of  theory, 
we  predict  them  (unsealed)  to  be  2284  cm-^,  1604  cm-l  and  1058  crn'T  Experimentally 
observed  isotopic  shifts  (C2H2  -  C2D2)  are  770  (2930-2160)  cm'^,  35  (1530-1495)  cm"^ 
and  365  (1375-1010)  cm‘^  respectively,  while  theoretically  predicted  (harmonic)  isotopic 
shifts  are  859  (3143-2284)  crn'K  33  (1637-1604)  cm-l  and  393  (1451-1058)  cm'l.  The 
agreement  between  experimental  and  theoretical  isotopic  shifts  is  very  satisfactory  and 
consistently  good  at  all  levels  of  theory  utihzed  in  this  paper. 

Transition  State  for  Triplet  Vinylidene  <->  Triplet  Acetylene  Isomerization 

The  Ci-symmetry  transition  state  for  isomerization  between  a  ^62  vinylidene  and 
b  ^Bu  trans-hent  acetylene  has  been  confirmed  by  vibrational  analysis  to  be  a  true  transition 
state  with  a  Hessian  index  of  one  at  all  levels  of  theory  employed  in  this  work  (see  Table 
III).  According  to  a  potential  energy  distribution  (FED)  analysis  of  the  normal  mode  with 
one  imaginary  frequency,  the  reaction  coordinate  consists  of  the  CCHmigrating  bend  and  the 
CHmigrating  Stretch.  With  the  TZ(2df,2pd)  basis,  the  imaginary  frequency  was  determined 
to  be  1731/  cm‘^  by  the  SCF  method  and  1500/  cm'*  by  the  CISD  method.  At  the  same 
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levels  of  theory,  the  lowest  real  harmonic  vibrational  frequency  was  determined  to  be  836 
cm-l  and  777  crn'l,  respectively.  The  large  magnitudes  of  the  imaginary  frequency  and  the 
lowest  real  harmonic  vibrational  frequency  suggest  that  the  nature  of  this  stationary  point 
will  not  change  even  if  more  sophisticated  correlated  methods  are  employed. 

The  harmonic  vibrational  frequency  for  the  CC  stretch  is  1637  cm'l  for  a  ^62 
vinylidene  and  1640  cm'l  for  b  ^Bu  acetylene  at  the  TZ(2df,2pd)  CISD  level  of  theory;  the 
CC  stretch  mode  for  the  transition  state  connecting  these  two  minima  has  a  frequency  of 
only  1355  cm'l.  This  red-shift  of  about  285  cm-l  is  due  to  the  lengthening  of  the  CC  bond 
in  the  transition  state  (bond  order  one  and  a  half)  relative  to  the  doubly-bonded  minima. 
Not  surprisingly,  the  vibrational  frequency  associated  with  the  CHmigrating  stretching  mode 
is  also  greatly  decreased  in  the  transition  state  structure  due  to  the  ~0.lA  increase  of  that 
bond  length.  The  average  of  the  two  TZ(2df,2pd)  CISD  CH  stretching  frequencies  of 
a  3B2  vinylidene  is  3190  cm'l,  and  for  b  ^Bu  acetylene  the  corresponding  average  is  3256 
cm'l.  In  the  transition  state  structure,  however,  the  harmonic  vibrational  frequency  for  the 
CHmigrating  mode  is  2316  cm-i,  which  is  red-shifted  by  -900  cm'l.  These  frequency 
shifts  with  changes  in  bond  lengths  are  expected  from  Badger's  Rules.46 

Transition  State  for  cis-trans  Triplet  Acetylene  Isomerization 

Vibrational  analysis  of  this  Cj-symmetry  stationary  point  has  confirmed  the 
structure  to  be  a  genuine  transition  state  with  only  one  imaginary  vibrational  frequency  at  all 
levels  of  theory  (see  Table  IV).  The  reaction  coordinate  is  mainly  the  CCHinverting  bend. 
In  contrast  with  the  transition  state  discussed  in  the  previous  section,  there  is  almost  no 
CHinverting  stretch  contribution  to  the  reaction  coordinate.  The  magnitude  of  the 
imaginary  vibrational  frequency  is  very  stable  with  respect  to  basis  set  improvements,  and 
it  only  decreases  by  about  150  cm'l  for  CISD  results  relative  to  SCF  results.  For  example, 
with  the  TZ(2df,2pd)  basis  set,  SCF  theory  predicts  a  magnitude  of  1 147i  cm'l,  while  the 
CISD  method  predicts  995i  cm'l.  The  stability  of  this  imaginary  frequency  and  the  lowest 
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real  vibrational  frequency,  the  a"  out-of-plane  bending  mode,  strongly  suggest  that  the 
nature  of  this  stationary  point  will  not  change  even  with  more  sophisticated  ab  initio 
methods. 

As  previously  mentioned,  this  transition  state  structure  has  a  CC  double  bond.  It  is 
expected  to  have  a  CC  stretch  vibrational  frequency  similar  to  those  of  a  ^62  and  b  ^Bu 
acetylene.  In  fact,  the  TZ(2df,2pd)  CISD  CC  stretch  frequency  of  1633  cm-1  for  this 
transition  state  structure  is  about  the  same  as  the  CC  stretch  frequencies  of  1660  cm'^  and 
1640  cm’l  predicted  for  a  ^B2  and  b  ^Bu  acetylene  with  an  equivalent  level  of  theory. 20 
Furthermore,  the  other  real  frequencies  have  magnitudes  similar  to  the  corresponding 
vibrational  frequencies  in  the  two  lowest  triplet  states  of  acetylene. 

Energetics 

d  ^82  Vinylidene 

In  comparison  with  results  published  by  Ervin  et  al.,^  the  ab  initio  theoretical 
predictions  for  the  (X  1  A]  -  a  3B2)  vinylidene  energy  separation  are  quite  reasonable.  At  all 
correlated  levels  these  results  (see  Table  VII)  fall  within  about  6  kcal/mole  of  the 
experimentally  reported  value.  The  theoretically  predicted  values  approach  the 
experimentally  reported  value  from  below  in  all  cases  as  the  basis  set  size  is  increased. 
Inclusion  of  electron  correlation,  also  tends  to  increase  the  predicted  singlet-triplet 
splittings.  As  improvements  in  the  treatment  of  electron  correlation  are  made  using  the 
same  basis  set,  the  predicted  value  of  the  splitting  seems  to  converge  toward  a  single  value. 
This  would  seem  to  indicate  that  the  CCSD(T)  single-point  energy  evaluations  at  the  CISD 
optimized  geometries  recover  most  of  the  electron  correlation  energy  for  both  of  these  states 
of  vinylidene.  Our  results  at  the  TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory  agree 
remarkably  well  with  experimentally  reported  values.  Specifically,  our  study  obtains  the 
result  Te(X  *Ai  -  a  3B2)  =  45.58  kcal/mole  compared  with  the  experimental  valued  of  46.8 
±  0.5  kcal/mole.  Likewise,  our  theoretical  prediction  of  To(X  ^Ai  -  a  ^82)  at  that  level  of 
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theory  is  46.56  kcal/mole,  while  the  experimentally  reported  value^  is  47.60  ±  0.14 
kcal/mole.  Experimental  and  theoretical  agreement  to  within  1  kcal/mole  is  generally 
considered  excellent.  This  instills  confidence  that  all  of  our  energetic  predictions  at  the 
TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory  should  be  within  1.5  kcal/mole  of  the  true 
values. 

As  shown  in  the  last  column  of  Table  VI,  the  "triplet-triplet"  splitting  between 
a  ^B2  vinylidene  and  a  ^62  acetylene  is  in  the  neighborhood  of  1-3  kcal/mole.  The 
classical  energy  difference  increases  slightly  with  basis  set  improvements  [e.g.  2.4 
kcal/mole  for  DZP  SCF  to  2.6  for  TZ(2df,2pd)  SCF]  and  generally  decreases  by  including 
electron  correlation  [e.g.  0.9  for  DZP  CCSD(T)//CISD].  The  best  value  predicted  in  this 
work,  with  ZPVE  corrections,  is  AEo(a  ^62  vinylidene  -  a  ^62  acetylene)  =  2.2  kcal/mole 
for  TZ(2df,2pd)  CCSD(T)//CISD. 

Triplet  Vinylidene  O  Triplet  Acetylene  Isomerization 

The  theoretically  predicted  energies  for  the  two  minima  involved  in  this 
isomerization  reaction,  determinable  from  Table  VI,  differ  by  5-6  kcal/mole.  With  ZPVE 
corrections,  at  the  TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory,  this  energy  difference  is 
6.1  (8.3  -  2.2)  kcal/mole.  The  reaction,  although  it  is  slightly  endothermic,  is  thus  nearly 
thermoneutral. 

With  the  DZP  basis,  the  classical  barrier  to  isomerization  via  1,2-hydrogen  shift  for 
the  SCF  method  is  60.7  (63.1  -  2.4)  kcal/mole,  but  it  drops  by  almost  4  kcal/mole  to  56.9 
(58.4  -  1.5)  kcal/mole  for  the  CISD  method  and  drops  further  to  53.2  (54.1  -  0.9) 
kcal/mole  for  the  CCSD(T)//CISD  method.  Basis  set  improvements  are  generally 
accompanied  by  decreases  in  the  predicted  value  of  the  classical  barrier.  For  example,  with 
the  CCSD(T)//CISD  method,  improvement  of  the  basis  set  quality  from  DZP  to 
TZ(2df,2pd)  leads  to  a  drop  in  the  classical  barrier  by  over  one  kcal/mole  to  52.1  (53.6  - 
1.5)  kcal/mole.  Finally,  for  all  basis  sets,  the  ZPVE  correction  lowers  the  barrier 
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significantly.  At  the  TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory,  with  CISD  ZPVE 
corrections,  the  barrier  to  isomerization  via  1,2-hydrogen  migration  is  47.9  (50.1  -  2.2) 
kcal/mole.  For  the  reverse  reaction,  the  barrier  is  41.8  (50.1  -  8.3)  kcal/mole. 

The  barrier  to  isomerization  predicted  for  the  forward  reaction  is  much  higher  than 
the  analogous  barrier  for  the  ground  state  vinylidene  reaction  of  1.6  kcal/mole.^  In  the  case 
of  the  singlet  ground  state,  facile  rearrangement  (admittedly,  with  extensive  rehybridization 
of  the  MO's)  can  take  place  by  migration  of  the  H  atom  toward  the  unoccupied  in-plane 
p-type  orbital  on  the  carbene  carbon.  On  the  other  hand,  in  the  case  of  a  ^62 
rearrangement,  the  p-type  MO  is  singly-occupied.  The  resultant  repulsive  interactions 
make  the  in-plane  migration  an  extremely  high-energy  pathway.  This  forces  the  triplet 
1,2-hydrogen  migration  to  take  place  via  another  pathway.  At  the  DZP  SCF  level  of 
theory,  the  Cs-symmetry-constrained  higher-order  saddle  point  lies  18.4  kcal/mole  above 
the  Cl -symmetry  transition  state. 

Transition  State  for  cis-trans  Triplet  Acetylene  Isomerization 

Based  upon  single-point  energies  evaluated  using  the  Cl  method  with  a  basis  set 
slightly  larger  in  size  than  our  DZP,  Demoulin'^^  suggested  that  the  classical  barrier  to 
cis-trans  triplet  acetylene  interconversion  by  internal  rotation  about  the  CC  bond  is  about 
23.1  kcal/mole  and  18.9  kcal/mole  for  the  reverse  reaction.  As  discussed  in  the  geometry 
section  above,  a  corresponding  stationary  point  along  the  internal  rotation  coordinate  was 
located  at  the  DZP  SCF  level  of  theory  with  a  total  energy  of  -76.671923  hartrees.  This 
places  it  35.3  kcal/mole  in  energy  above  a  ^32  acetylene,  although  it  is  not  a  true  transition 
state.  By  comparison,  the  classical  barrier  to  interconversion  along  an  inversion  pathway  is 
only  17.3  kcal/mole  at  the  DZP  SCF  level  of  theory.  Inversion  is  thus  the  energetically 
preferred  process  for  the  cis-trans  isomerization. 

As  shown  in  the  third  column  of  Table  VI,  the  barrier  to  cis-trans  isomerization  via 
in-plane  inversion  decreases  with  both  basis  set  improvements  and  inclusion  of  electron 
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correlation  effects.  According  to  arguments  similar  to  those  presented  in  the  previous  two 
sections,  the  predicted  barriers  should  converge  to  within  a  couple  kcal/mole  of  the  true 
value.  At  the  TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory,  we  found  the  ZPVE  corrected 
barrier  to  be  13.0  kcal/mole. 

The  exothermic  reverse  reaction,  trans-cis  triplet  acetylene  interconversion  by 
inversion,  has  a  markedly  lower  barrier  to  rearrangement.  Our  best  prediction  of  the  trans- 
cis  triplet  acetylene  isomerization  barrier  is  4.7  (13.0  -  8.3)  kcal/mole  with  ZPVE 
corrections  at  the  TZ(2df,2pd)  CCSD(T)//CISD  level  of  theory.  Although  this  barrier  is 
quite  small,  it  is  unlikely  to  disappear  entirely  at  higher  levels  of  theory. 

Conclusion 

High-level  ab  initio  molecular  electronic  structure  theory  has  been  applied  to 
investigate  the  a  ^62  electronic  state  of  vinylidene,  the  transition  state  for  isomerization  of  a 
3B2  vinylidene  to  b  ^Bu  trans-bent  acetylene,  and  the  transition  state  for  interconversion 
between  a  ^B2  cw-bent  and  b  ^Bu  trans-bent  acetylene. 

The  a  3B2  state  of  vinylidene  is  confirmed  to  be  a  doubly-bonded  minimum-energy 
structure  at  all  levels  of  theory.  Predicted  harmonic  vibrational  frequencies  agree  well  with 
available  experimental  fundamental  frequencies  for  a  ^B2  vinylidene  and  for  the  mono-  and 
di-  deuterated  isotopomers  of  a  3B2  vinylidene.  Theoretical  predictions  of  To(X  lAi  - 
a  ^B2)  for  vinylidene  are  also  in  excellent  agreement  with  the  results  of  negative  ion 
photoelectron  spectroscopy  experiments. 

A  nonplanar  bridged  structure  was  found  to  be  the  transition  state  for  the 
rearrangement  of  a  ^B2  vinylidene  to  b  ^Bu  trans-bent  acetylene.  The  barrier  to 
1,2-hydrogen  migration  was  found  to  be  47.9  kcal/mole  at  the  highest  level  of  theory, 
TZ(2df,2pd)  CCSD(T)//CISD. 

The  unimolecular  reaction  which  interconverts  a  ^B2  cw-bent  and  b  ^Bu  trans-bent 
acetylene  was  found  to  occur  via  a  planar  inversion  pathway  rather  than  by  rotation  about 
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the  CC  bond  axis.  The  barrier  for  this  inversion  process  was  shown  to  be  only  13.0 
kcal/mole  at  the  highest  level  of  theory. 
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Appendix  A 

The  TZ(2df,2pd)  CISD  harmonic  vibrational  frequencies  for  the  3B2  state  of 
acetylene  are  3196  cm'l,  1660  cm'*,  and  849  cm'l  for  the  ai  modes.  The  a2  mode 
frequency  is  839  cm'l’  and  the  b2  modes  predictions  are  3160  cm'l  and  1 133  cm‘l.  Those 
for  the  ^Bu  state  of  acetylene  are  3265  cm‘l,  1640  cm'f ,  and  1 124  cm‘1  for  the  ag  modes, 
937  cm'l  for  the  au  mode,  and  3247  cm'*  and  791  cm'*  for  the  bu  modes. 

For  the  DZP  and  TZ2P  basis  sets,  the  CCSD//CISD  energies  of  X  *Ai  vinylidene 
are  -77.038632  and  -77.078175,  respectively.  For  CCSD(T)//CISD  with  those  basis  sets, 
the  energies  are  -77.047967  and  -77.090869.  With  the  TZ(2df,2pd)  basis  set,  the  SCF 
energy  at  the  optimized  geometry  is  -76.795671.  Harmonic  vibrational  frequencies  at  this 
level  of  theory  are  3363  cm'*,  3274  cm'*,  1819  cm'*,  1352  cm'*,  912  cm'*  and  523  cm'*. 
The  TZ(2df,2pd)  CISD  energy  is  -77.081828  hartrees.  Harmonic  vibrational  frequencies 
for  this  level  of  theory  are  3294  cm'*,  3201  cm'*,  1732  cm'*,  1265  cm'*,  798  cm'*  and 
376  cm-*.  The  TZ(2df,2pd)  CISD+Q//CISD,  CCSD//CISD  and  CCSD(T)//CISD  energies 
are  -77.110428  hartrees,  -77.105181  hartrees  and  -77.119439  hartrees,  respectively. 
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Figure  1:  Predicted  Structure  of  a  ^62  Vinylidene. 
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Figure  2:  Predicted  Structure  of  Ci -symmetry  Transition  State  for  the  %2  Vinylidene 
3Bu  trans-hsai  Acetylene  Isomerization  Reaction. 


Figure  3:  Predicted  Torsional  Angle  of  the  Ci-symmetry  Transition  State  for  the  ^62 
Vinylidene  o  ^Bu  trans-h&cA  Acetylene  Isomerization  Reaction. 
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DZPSCF  1.3233 A 


TZ2PSCF  1.3120  I.O6OIA 

TZ(2df,2pd)  SCF  1.3114  1  0529 

1.0534 

DZPCISD  1.3442 

TZ2PCISD  1.3281  1.0696 

TZ(2df,2pd)  CISD  1.3248  1.0587 
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Figure  4:  Predicted  Structure  of  Cs-symmetry  Transition  State  for  the  ^Bu  trans-bent 
3B2  c/s-bent  Acetylene  Isomerization  Reaction. 


Table  I.  Theoretical  predictions  of  the  total  energy  (in  hartrees),  dipole  moment  (in  debye),  harmonic  vibrational  frequencies  (in  cm'l), 
infrared  (IR)  intensities  (in  km  mol'^)  and  zero  point  vibrational  energy  (ZPVE)  (in  kcal  mol"^)  for  the  a  ^62  state  of  vinylidene. 
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Table  III.  Theoretical  predictions  of  the  total  energy  (in  hartrees),  dipole  moment  (in  debye),  harmonic  vibrational  frequencies  (in 
cm-l),  and  zero  point  vibrational  energy  (ZPVE)  (in  kcal/mole)  for  the  Ci-symmetry  transition  state  between  a  ^62  vinylidene  and  b  ^Bu 
acetylene. 
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Table  IV.  Theoretical  predictions  of  the  total  energy  (in  hartrees),  dipole  moment  (in  debye),  harmonic  vibrational  frequencies  (in 
cm-1),  and  zero  point  vibrational  energy  (ZPVE)  (in  kcal/mole)  for  the  Cg-symmetry  transition  state  between  a  ^62  acetylene  and  b  ^Bu 
acetylene. 
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CISD+Q,  CCSD  and  CCSD(T)  energies  were  evaluated  at  the  CISD  optimized  geometries  for  the  same  basis  set. 


Table  VI.  Relative  Energies^  (in  keal/mole)  for  the  triplet  acetylene/vinylidene  potential  energy  hypersurface. 
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TZ2PCCSD(T)  0.00  14.48  (13.10)  7.59  (7.91)  54.56  (50.95)  0.91  (1.70) 

TZ(2df,2pd)  CCSD(T)c _ OOO _ 14.33  (13.04)  8.03  (8.27)  53.62  (50.10) _ 1.48  (2.22) 


^Numbers  in  parentheses  include  zero  point  vibrational  energy  (ZPVE)  corrections.  CISD+Q,  CCSD  and  CCSD(T)  energies  were 
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Table  VII.  Theoretical  predictions  for  the  X  ^Ai-  a  3B2  vinylidene  energy  separation  (in 
kcal/mole)®’*’. 


Level  of  Theory 

Te.(XiAi-a3B2) 

To(X  lAi-a3B2) 

DZPSCF 

31.37 

32.31 

TZ2P  SCF 

32.79 

33.76 

TZ(2df,2pd)  SCF 

32.89 

33.86 

DZP  CISD 

40.70 

41.56 

TZ2P  CISD 

42.36 

43.30 

TZ(2df,2pd)  CISD 

43.75 

44.73 

DZP  CISD+Q 

41.45 

42.31 

TZ2P  CISD+Q 

43.53 

44.47 

TZ(2df,2pd)  CISD+Q 

45.16 

46.14 

DZP  CCSD 

41.21 

42.06 

TZ2P  CCSD 

43.18 

44.12 

TZ(2df,2pd)  CCSD 

44.76 

45.74 

DZP  CCSD(T) 

41.76 

42.61 

TZ2P  CCSD(T) 

43.93 

44.87 

TZ(2df,2pd)  CCSD(T) 

45.58 

46.56 

Expt.*’ 

46.8  ±  0.5 

47.60  ±  0.14 

aCISD+Q,  CCSD  and  CCSD(T)  energies  were  determined  at  the  CISD  optimized 


geometries,  and  the  CISD  ZPVE  corrections  with  the  same  basis  set  were  used.  Where 
possible,  energies  for  X  ^ A]  vinylidene  were  taken  from  reference  7.  The  rest  were 
obtained  in  this  work  (see  Appendix  A). 

•’The  experimentally  measured  values  are  taken  from  reference  6. 


Chapter  III 


The  X  AlOH  -  X  HAIO  Isomerization  Potential  Energy  Hypersurface’ 


’  G.  Vacek,  B.  J.  DeLeeuw  and  H.  F.  Schaefer  7.  Chem.  Phys.  98,  8704  (1993). 
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Ab  initio  molecular  electronic  structure  theory  has  been  used  to 
study  the  AlOH-HAlO  unimolecular  isomerization  reaction  on  the  singlet 
ground  state  potential  energy  surface.  Electron  correlation  effects  were 
included  via  configuration  interaction  and  coupled-cluster  methods.  Basis 
sets  as  complete  as  triple  zeta  plus  two  sets  of  polarization  functions  and  a 
set  of  higher  angular  momenrnm  functions  lTZ(2df,2pd)]  were  employed. 
The  classical  barrier  for  hydrogen  migration  from  X  HAIO  to  X  AlOH  is 
predicted  to  be  38.4  kcal  mol'*  using  the  TZ(2df,2pd)  basis  set  with  the 
coupled-cluster  method  including  all  single  and  double  excitations  with  the 
effect  of  connected  triple  excitations  included  perturbatively  [CCSD(T)]. 
After  correction  for  zero-point  vibrational  energies  (ZPVEs),  an  activation 
energy  of  36.6  kcal  mok*  is  obtained.  The  AE  for  isomerization  is  -42.2 
(-40.5  with  ZPVE  correction)  kcal  mok*  at  the  same  level  of  theory.  The 
dipole  moments  of  HAIO  and  AlOH  in  their  equilibrium  geometries  are 
4.525  and  1.040  debye,  respectively,  at  the  same  level  of  theory.  A 
comparison  is  also  made  between  a  theoretically  predicted  harmonic 
vibrational  frequency  and  a  recently  determined  experimental  fundamental 

frequency  for  X  AlOH. 


Introduction 

Metal  hydroxides  are  often  the  result  of  metal  atom  insertion  or  elimination  reactions 
with  water  in  environments  undergoing  metal  ablation,  sputtering  or  combustion.  The 
simplest  of  these,  metal  monohydroxides,  form  a  class  of  molecules  upon  which  to  test 
elementary  concepts  of  chemical  bonding.  Nonlinearity  of  metal  hydroxides  is  thought  to 
indicate  the  extent  to  which  the  metal-hydroxide  bond  is  covalent  rather  than  ionic. 
Covalent  bonding  tends  to  favor  bent  structures  (e.g.  H2O,  BOH*),  while  ionic  bonding 
favors  linear  structures  (e.g.  CaOH,^  SrOH,^  BaOH^).  MgOH,^  which  is  quasilinear, 
may  be  considered  to  have  an  ionic  bond  with  covalent  features.  Evidence  suggests  that  the 
linearity,  nonlinearity,  or  quasilinearity  of  a  molecule  may  be  qualitatively  predicted  from 
simple  electronegativity  arguments.  Atoms  such  as  B  and  Cu  with  electronegativities  near 
that  of  hydrogen  (2.2)  form  bent  monohydroxides. *>^  Atoms  with  electronegativities  near 
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or  below  1,  like  most  alkali  metals^.^  and  alkaline  earths, 2-4  turn  out  to  be  linear.  Mg  has 
an  electronegativity  between  these  two  extremes  (1.31),  and  MgOH  turns  out  to  be 
quasilinear.5  a1  has  an  electronegativity  of  1.61,  so  AlOH  might  also  expected  to  be 
quasilinear. 

Recently,  the  first  observation  of  an  electronic  spectrum  for  X  AlOH  has  been 
reported  by  Pilgrim,  Robbins  and  Duncan.^  The  presumed  structure  and  bonding  of  the 
system  was  also  discussed.  Since  they  produced  the  subject  of  their  research  through 
aluminum  insertion  into  water,  it  was  initially  unclear  whether  their  spectrum  should  be 
assigned  to  the  AlOH  or  HAIO  isomer.  Previous  theoretical  studies,>o,i i 
mass-spectrometric  analyses,i2  and  infrared^^  (IR)  and  optical*^  assignments  did  not 
provide  sufficient  information  about  the  relative  stabihties  and  the  vibrational  frequencies  of 
these  two  species  for  them  to  make  a  satisfactory  assignment.  Here  we  present  results  of 
our  inquiry  into  the  isomerization  reaction  from  X  AlOH  to  X  HAIO  and  the  associated 
potential  energy  surface  (PES).  Some  of  the  data  reported  herein  was  crucial  to  Pilgrim 
et  ai's^  assignments.  We  are  confident  that  it  will  be  helpful  to  future  experimental  studies 
of  these  two  interesting  isomers. 

Theoretical  Methods 

Three  different  basis  sets  were  used  in  this  research.  The  double  zeta  plus 
polarization  basis  (DZP)  is  a  standard  Huzinaga-Dunning-Hayi5,i6  double-^  set  of 
contracted  Gaussian  functions  augmented  by  a  set  of  six  Cartesian  d-like  polarization 
functions  on  oxygen  [ad(0)  =  0.85]  and  aluminum  [ad(Al)  =  0.40]  and  a  set  of  p-type 
polarization  functions  on  hydrogen  [ap(H)  =  0.75].  The  contraction  scheme  for  the  DZP 
basis  is  thus: 

Al(lls7pld/6s4pld);  0(9s5pld/4s2pld);  H  (4slp/2slp). 

The  triple  zeta  plus  double  polarization  basis  (TZ2P)  is  triple-^  in  the  valence  region  and 
adds  two  sets  of  polarization  functions  per  atom  [ad(Al)  =  0.8, 0.2;  ad(0)  =  1.7,  0.425; 
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ap(H)  =  1.5,  0.375  ].  For  hydrogen  and  oxygen,  standard  Huzinaga-Dunningi^.n  sets  of 
contracted  Gaussian  functions  were  used.  For  aluminum,  McLean  and  Chandler's** 
contraction  of  Huzinaga's  12s9p  primitive  Gaussian  set  was  used.  The  contraction  scheme 

for  the  TZ2P  basis  is  as  follows: 

A1  (12s9p2d/6s5p2d);  O  (10s6p2d/5s3p2d);  H  (5s2p/3s2p). 

The  TZ(2df,2pd)  basis  was  formed  by  adding  a  set  of  higher  angular  momentum 
polarization  functions  [aKAl)  =  0.25;  tXf(0)  =  1.4;  and  cXd(H)  =  1.0]  to  the  TZ2P  basis. 

The  restricted  Hartree-Fock  self-consistent  field  (SCF)  method  was  initially  employed 
in  this  study.  The  effects  of  electron  correlation  were  taken  into  account  using  the  method 
of  configuration  interaction  including  all  single  and  double  excitations  from  an  SCF 
reference  wavefunction  (CISD),  the  coupled-cluster  method  including  all  single  and  double 
excitations  (CCSD),  and  the  CCSD  method  with  the  effect  of  connected  triple  excitations 
included  perturbatively  [CCSD(T)].  The  two  lowest  (aluminum  and  oxygen  Is-like)  SCF 
molecular  orbitals  were  constrained  to  be  doubly-occupied  and  the  corresponding  two 
highest  virtual  orbitals  were  deleted  from  the  correlation  procedure.  With  the  TZ(2df,2pd) 
basis,  this  resulted  in  a  total  of  72,311  configurations  in  C2v  symmetry  and  138,919 
configurations  in  Cs  symmetry. 

The  structures  of  ground  state  AlOH,  HAIO,  and  the  transition  state  connecting  them 
via  hydrogen  migration  were  fully  optimized  using  closed-shell  analytic  gradient  techniques 
at  the  SCF,i9-2i  ciSD,22-24  CCSD,25  and  CCSD(T)26  levels  of  theory.  Harmonic 
vibrational  firequencies  were  obtained  using  analytic  SCF  second-derivative  techniques^^'^^ 
and  central  finite  differences  of  analytic  CISD,  CCSD,  and  CCSD(T)  gradients. 

Relative  energies  designated  CISD+Q  were  obtained  for  the  various  isomers  studied 
here  by  adding  the  Davidson  corrections**’^*  for  unlinked  quadruple  excitations  to  the  CISD 


energy. 
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Discussion 

Geometry 

Walsh's  rules32  predict  that  both  AlOH  and  HAIO  should  have  linear  ground  state 
structures,  just  like  the  isovalent  molecules  HCN  and  HNC.  Furthermore,  it  can  be 
expected  from  the  same  rules  that  the  lowest  excited  states  of  AlOH  and  HAIO  should  be 
bent  and  have  similar  bond  angles.  However,  even  in  the  simple  case  of  the  excited  states 
of  HCN  and  HNC,  it  has  been  shown^s  that  Walsh  concepts  are  not  well  followed. 

As  shown  in  Figure  1,  the  use  of  a  DZP  quality  basis  set  leads  to  the  prediction  of  a 
linear  X  AlOH  structure  with  all  theoretical  methods  employed  here.  Previous  theoretical 
resultsiO’ii  also  predicted  AlOH  to  be  linear  using  basis  sets  similar  in  extent  to  our  DZP 
basis.  However,  with  the  more  complete  TZ2P  basis  set,  we  predict  a  geometry  that  is 
quite  bent  (-151°).  This  qualitative  difference  between  the  DZP  and  TZ2P  equilibrium 
X  AlOH  geometries  implies  that  the  use  of  still  larger  basis  sets  is  desirable.  It  is  not 
uncommon  for  a  basis  set  with  extensive  character  in  the  spd  space  but  no  higher  angular 
momentum  functions  to  prefer,  somewhat  erroneously,  a  bent  geometry  over  a  linear 
one. 34.35  In  such  cases,  further  study  with  basis  sets  which  do  incorporate  /-type 
functions  will  usually  lead  to  an  accurate  prediction  of  both  geometry  and  harmonic 

vibrational  frequencies.36.37 

Thus,  our  best  prediction  of  the  equilibrium  geometry  of  AlOH  is  at  the  TZ(2df,2pd) 
CCSD(T)  level  of  theory,  where  the  bond  angle  is  predicted  to  be  162.6°.  At  this  level  of 

o 

theory,  we  see  essentially  a  hydroxide  group  with  a  tight  single  bond  (re(OH)  =  0.951  A). 
The  experimentally  estimated  value  of  re(OH~)  lies  within  0.001  A38  of  the  re(OH)  value, 
0.970  A.39  This  "hydroxide  group"  is  ionically  bonded  to  the  partially  positive  A1  atom 
(re(AlO)  =  1.682  A).  Plots  of  the  9a'  molecular  orbital  (MO)  show  that  it  is  essentially  a 
lone  pair  of  electrons  in  an  A1  3^  orbital  slightly  polarized  along  the  molecular  axis  by  the 
A1  3pz  orbital  to  produce  a  weak  sigma  type  interaction  between  this  "lone  pair"  and  the 
hydroxide  group. 
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It  is  important  to  note  that  despite  the  bent  equilibrium  geometry  of  AlOH,  it  can  at 
best  be  referred  to  as  quasilinear.  For  many  of  the  levels  of  theory  utilized  in  this  paper, 
barriers  to  linearity  were  calculated.  In  all  cases,  the  predicted  barrier  to  linearity  was  less 
than  1  kcal  mol'K  Thus  the  ZPVE  is  greater  than  the  linearity  barrier.  In  many 
experimental  set-ups,  the  "averaged"  linear  geometry  would  be  observed.  The  predicted 
bond  angle  is  significantly  larger  than  that  predicted  for  X  BOH  (0e  =  125  ),^  which  has  a 
truly  bent  structure.  This  stmctural  difference  between  AlOH  and  BOH  is  in  agreement 
with  the  qualitative  predictions  based  upon  electronegativities  mentioned  in  the 
introduction. 

The  strucmre  of  X  HAIO,  depicted  in  Figure  2,  was  predicted  to  be  linear  at  all  levels 
of  theory,  in  accordance  with  Walsh's  rules.32  The  best  results,  at  the  TZ(2df,2pd) 
CCSD(T)  level  of  theory,  show  re(AlH)  =  1.562  A  and  re(AlO)  =  1.608  A.  The  AlH  bond 
distance  is  somewhat  short  compared  to  the  1.6478  A  bond  length  of  diatomic  AlH.^^ 
Presumably,  this  is  due  to  the  greater  s  character  within  the  AlH  bond  of  HAIO  as 
compared  to  that  of  the  diatomic.  According  to  Huber  and  Herzberg,^^  the  bond  length  of 
the  AlO  diatomic  molecule  (with  a  bond  order  of  2.5)  is  1.6179  A. 

Figure  3  shows  the  transition  state  for  the  hydrogen  migration  from  AlOH  to  HAIO. 
The  hydrogen  is  closer  to  A1  than  O,  which  agrees  nicely  with  Hammond’s  Postulate  for 
endothermic  reactions,^®  since  HAIO  is  much  higher  in  energy  than  AlOH.  The  best 
description  of  the  transition  state  structure  should  again  be  the  TZ(2df,2pd)  CCSD(T) 
results. 


Harmonic  Vibrational  Frequencies 

Typically,  the  harmonic  vibrational  frequencies  of  a  molecule  are  only  predicted 
accurately  at  levels  of  theory  where  the  predicted  geometry  is  qualitatively  correct.  We 
suspect  that  only  the  frequencies  predicted  with  the  TZ(2df,2pd)  basis  set  and  reasonable 
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levels  of  electron  correlation  will  thus  be  accurate  to  an  acceptable  degree.  Specifically,  we 
expect  the  TZ(2df,2pd)  CCSD(T)  level  of  theory  to  give  excellent  results.37 

At  the  TZ(2df,2pd)  CCSD(T)  level  of  theory,  our  predictions  of  the  harmonic 
vibrational  frequencies  of  AlOH,  found  in  Table  I,  are  4017, 845  and  155  cm-i  for  the  OH 
stretching,  AlO  stretching,  and  the  bending  frequencies.  The  only  experimentally  assigned 
fundamental  vibrational  frequency  in  the  gas  phase  is  that  of  Pilgrim.^  A  AG1/2  value  of 
895  cm-i  was  tentatively  assigned  to  the  AlO  stretching  mode.  This  is  higher  than  our 
prediction,  but  the  agreement  is  still  acceptable.  It  is  also  conceivable  that  Pilgrim  observed 
a  combination  band  of  the  AlO  stretch  and  the  AlOH  bend.  The  sum  of  our  harmonic 
vibrational  predictions  is  1000  cm'*,  nicely  above  Pilgrim's  895  cm'*  fundamental.  Our 
predictions  for  the  two  harmonic  stretching  vibrational  frequencies  are  also  in  satisfactory 
agreement  with  the  fundamental  stretching  frequencies  of  Hauge  et  al.*^  (Vi  =  3790,  V2  = 
810  cm'*),  which  were  measured  in  a  noble  gas  matrix.  The  disagreement  between  these 
two  experimental  V2  assignments  is  somewhat  surprising. 

For  the  HAIO  conformer.  Table  E,  the  TZ(2df,2pd)  CCSD(T)  harmonic  vibrational 
frequencies  are  predicted  to  be  2002, 1056  and  400  cm'*  for  the  AlH  stretching  mode,  the 
AlO  stretching  mode  and  the  bending  mode,  respectively.  These  predicted  values  do  not 
match  well  with  Pilgrim’s^  results,  thus  supporting  his  assignment  to  AlOH  rather  than 
HAIO. 

There  are  no  experimental  vibrational  frequency  assignments  for  HAIO.  Huber  and 
Herzberg39  jjgt  experimental  harmonic  vibrational  frequency  assignments  for  the  diatomic 
molecules  AlH  and  AlO  as  1683  and  979  cm'*,  respectively.  The  HAIO  molecule  has  an 
AlO  stretch  harmonic  vibrational  frequency  that  looks  just  slightly  perturbed  from  that  of 
diatomic  AlO. 

Harmonic  vibrational  frequencies  for  the  isotopically-substituted  ground  states  are 
presented  in  Tables  III  and  IV.  Comparison  between  predicted  harmonic  and 
experimentally  assigned  fundamentals  are  favorable  for  AlOD.  The  only  isotopic  shift  for 
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AlOD  that  is  determinable  from  experimental  results  is  AV2  =  15.1  cm-*  by  Hauge.*^  The 
predicted  isotopic  shift  agrees  well  at  all  levels  of  theory,  e.g.  18  cm-*  for  TZ(2df,2pd) 
CCSD(T).  Hauge*3  also  reports  that  V2(A1*80)  =  785.2  cm-*,  shifted  25.1  cm-*.  Our 
theoretical  predictions,  not  tabulated  here,  agree  well  with  this  experimental  assignment. 
For  TZ(2df,2pd)  CCSD(T)  we  predict  a)e(Al*^0)  =  817  cm-*,  a  shift  of  28  cm-*.  There  are 
no  experimental  isotopic  vibrational  assignments  for  DAIO,  but  we  are  confident  that  our 
predictions  are  good. 

Table  V  shows  the  harmonic  vibrational  frequencies  of  the  transition  state  for 
hydrogen  migration.  At  the  TZ(2df,2pd)  CCSD(T)  level  of  theory,  the  real  vibrational 
modes  have  the  magnitudes  1755  cm-*  and  909  cm-*,  while  the  imaginary  mode  has  a 
magnitude  of  9451  cm-*. 

Energetics 

Relative  energies  of  X  AlOH,  X  HAIO,  and  the  transition  state  are  given  in  Table  VI. 
At  the  highest  level  of  theory  used  here  [TZ(2df,2pd)  CCSD(T)],  AlOH  is  more  stable  than 
HAIO  by  42.2  kcal  mob*.  At  the  same  level  of  theory  the  classical  barrier  to  hydrogen 
migration  from  X  AlOH  to  X  HAIO  is  found  to  be  80.5  kcal  mob*,  while  zero-point 
vibrational  energy  (ZPVE)  corrections  give  an  activation  energy  of  77.2  kcal  mob*.  Thus, 
AlOH  is  extremely  stable  with  respect  to  isomerization  once  produced.  Although  the  HAIO 
to  AlOH  isomerization  reaction  is  strongly  exothermic  with  AE  =  -42.2  (-40.5  with  ZPVE) 
kcal  mob*,  the  barrier  height  for  X  HAIO  to  X  AlOH  isomerization  is  38.4  (36.6  with 
ZPVE)  kcal  mob*.  If  it  were  generated,  HAIO  should  be  observable  because  of  its  kinetic 
stabihty  with  respect  to  isomerization. 

Interestingly,  the  PES  in  the  related  case  of  X  HBO  o  X  BOH*’**  isomerization  is 
almost  the  exact  opposite  of  that  predicted  here  for  HAIO  and  AlOH.  That  is,  the  trivalent 
HBO  isomer  is  ~40  kcal  mob*  more  stable  than  the  monovalent  BOH  isomer,  and  the 
transition  state  connecting  them  is  ~50  kcal  mob*  above  BOH.*  In  many  ways,  first-row 
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elements  such  as  boron  are  oddities  and  stand  apart  from  the  heavier  elements  in  their 
respective  groups.  When  moving  down  the  periodic  table  within  Group  13  (B,  Al,  Ga,  In, 
Tl),  the  monovalent  state  becomes  progressively  more  stable.^^  One  explanation  for  this 
trend  and  the  predicted  monovalent  AlOH  ground  state  is  the  inert  5-pair  effect. 
Kutzelnigg42  has  noted  that  the  2s  and  2p  orbitals  of  first-row  elements  are  uniquely  well 
suited  for  hybridization  because  they  have  similar  spatial  extent.  First-row  elements  have 
only  5-type  core  electrons.  As  a  result,  their  valence  p  orbitals  are  relatively  compact. 
Elements  in  higher  rows  have  cores  that  contain  both  s-  and  p-type  electrons;  their  valence 
p  orbitals  are  substantially  more  diffuse  than  their  valence  s  orbitals,  and  hybridization  is 
much  less  efficient. 

Recall  that  the  chemistry  of  lower  Group  14  elements  (Si,  Ge,  etc.)  is  often  radically 
different  from  that  of  carbon;  this  has  been  explained  in  terms  of  increasing  stability  of  the 
divalent  state  while  moving  downward  within  the  group.  The  divalent  state  stabilization 
energy  (DSSE),  defined  to  be  the  difference  between  the  first  and  second  bond  dissociation 
energies,  is  a  convenient  way  to  express  this  trend.  For  example,  DSSE(CH2)  =  -5.6 
kcal  mol-i,  while  DSSE(SiH2)  =  +19.3  kcal  mofi  and  DSSE(GeH2)  =  +25.8  kcal  mokf 
The  usefulness  of  DSSEs  with  respect  to  Group  14  chemistry  has  been  discussed 
elsewhere  by  Grev.'^^  Using  G2  theory  thermochemical  data,'^^  the  equivalent 
"monovalent  state  stabilization  energy"  may  be  determined  for  both  BH  (-28.8  kcal  mok') 
and  AlH  (+39.6  kcal  mol-')-  Although  this  treatment  of  the  thermochemistry  is  hardly 
complete,  it  does  provide  a  strong  indication  that  the  monovalent  structure  (MOH)  should 
be  much  more  favorable  relative  to  the  trivalent  structure  (OMH)  in  the  case  of  Al  as 
opposed  to  B. 

Conclusion 

The  singlet  ground  state  PES  for  the  AlOH-HAlO  isomerization  reaction  has  been 
studied  with  ab  initio  molecular  electronic  structure  theory.  AlOH  was  found  to  be 
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quasilinear  with  an  equilibrium  bond  angle  of  162.6°  contrary  to  Walsh's  rules. 
Predictions  herein  for  AlOH  compare  favorably  with  recently  determined  experimental 
results.  Agreement  is  sufficient  to  confirm  the  experimental  assigmnents.  HAIO  was  linear 
at  all  levels  of  theory  in  agreement  with  Walsh. 

At  the  [TZ(2df,2pd)]  CCSD(T)  level  of  theory,  the  classical  barrier  for  X  HAIO  to 
X  AlOH  rearrangement  is  predicted  to  be  38.4  kcal  mob^  while  it  is  36.6  kcal  mob’  with 
ZPVE  corrections.  The  AE  for  this  isomerization  reaction  is  -42.2  (-40.5  with  ZPVE)  kcal 
mol'i  at  that  level  of  theory.  The  AlOH  conformation  lies  lowest  on  the  AlOH-HAlO  PES, 
in  stark  contrast  with  the  HBO  global  minimum  on  the  BOH-HBO  PES.  This  apparent 
lack  of  periodicity  within  Group  13  can  be  explained  in  terms  of  the  increasing  stability  of 
the  monovalent  state  for  lower  elements  in  that  group. 
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Figure  1.  Predicted  equilibrium  geometry  for  X  AlOH. 
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Figure  2.  Predicted  equilibrium  geometry  for  X  HAIO. 
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Figure  3.  Predicted  structure  for  Cs-symmetry  transition  state  for  the  X  AlOH 


X  HAIO  isomerization  reaction. 
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Vi  =  3790  V2  =  810.3 


DZP  equilibrium  geometries  are  linear,  while  TZ2P  and  TZ(2df,2pd)  are  bent. 


57 


o 

> 

o 

(L> 

Oh 

ua 

8 

c/T 

o 

H 

% 

XJ 

VO 


00 

00 


cn 

I 

§ 

OO*' 

m 

rt 

O 

00 

cn 


cs 

VO 


r- 


cn 

j 

B 

a 

D 

(D 

a 

W 

a 

+ 

Q 

C4) 

t— ( 

U 


C3 


Table  II.  Theoretical  predictions  of  the  total  energy  (in  hartrees),  dipole  moment  (in  debye),  harmonic  vibrational  frequencies  (in  cm‘1), 
infrared  (IR)  intensities  (in  km  mol'i)  and  zero  point  vibrational  energy  (ZPVE)  (in  kcal  mol'*)  for  X  HAIO. 
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CISD+Q  Energies  are  -317.701  586,  -317.779  378,  and  -317.810418  hartrees,  respectively. 
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Table  III.  Theoretical  predictions  of  the  harmonic  vibrational  frequencies  (in  cm*!)  for 
AlOD,  the  deuterated  isotopomer  of  ground  state  of  AlOH. 

Harmonic  Vibrational  Frequency 


Level  of  Theory  C0i(ai/a')  a)2(ai/a')  C03(ei/a’) 


DZP  SCF 

3162 

880 

113 

TZ2P  SCF 

3134 

853 

194 

TZ(2df,2pd)  SCF 

3147 

866 

117 

DZP  CISD 

3036 

864 

101 

TZ2PCISD 

3023 

832 

149 

TZ(2df,2pd)  CISD 

3034 

850 

91 

DZPCCSD 

2980 

855 

77 

TZ2P  CCSD 

2955 

817 

167 

TZ(2df,2pd)  CCSD 

2967 

836 

_  115 

DZP  CCSD(T) 

2962 

851 

60 

TZ2PCCSD(T) 

2922 

808 

166 

TZ(2df,2pd)  CCSD(T) 

2932 

827 

118 

Expt.7 

Expt.5 

V  =  2600  ±  100 

V  =  750  ±  100 

V  =  795.2 

V  =  240  ±  50 
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Table  IV.  Theoretical  predictions  of  the  harmonic  vibrational  frequencies  (in  cm-1)  for 
DAIO,  the  deuterated  isotopomer  of  ground  state  of  HAIO. 


Harmonic  Vibrational  Frequency 


Level  of  Theory 

COi(ai) 

(02(ai) 

C03(ei) 

DZP  SCF 

1566 

1186 

377 

TZ2P  SCF 

1543 

1174 

370 

TZ(2df,2pd)  SCF 

1540 

1172 

372 

DZP  CISD 

1552 

1118 

347 

TZ2P  CISD 

1514 

1106 

341 

TZ(2df,2pd)  CISD 

1517 

1113 

343 

DZPCCSD 

1521 

1085 

328 

TZ2P  CCSD 

1475 

1070 

321 

TZ(2df,2pd)  CCSD 

1476 

1077 

_  325 

DZP  CCSD(T) 

1503 

1046 

312 

TZ2P  CCSD(T) 

1456 

1030 

303 

TZ(2df,2pd)  CCSD(T) 

1458 

1038 

306 

Table  V.  Theoretical  predictions  of  the  total  energy  (in  hartrees),  dipole  moment  (in  debye),  harmonic  vibrational  frequencies  (in  cm"'), 
and  zero  point  vibrational  energy  (ZPVE)  (in  kcal/mole)  for  the  Cs-symmetry  transition  state  connecting  X  AlOH  and  X  HAIO. 
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CISD+Q  Energies  are  -317.636  829,  -317.713  007,  and  -317.742  528  hartrees,  respectively. 
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Table  VI.  Relative  Energies^  (in  kcal/mole)  for  the  ground  HAIO-AIOH  potential  energy 
hypersurface. 


Level  of  Theory 

X  AlOH 

Cs-symmetry 

transition  state 

X  HAIO 

Barrier 

Height 

DZP  SCF 

0.00 

112.33  (108.64) 

54.72  (53.02) 

57.60  (55.61) 

TZ2P  SCF 

0.00 

110.88  (107.33) 

52.63  (50.98) 

58.26  (56.35) 

TZ(2df,2pd)  SCF 

0.00 

111.71  (108.27) 

53.33  (51.79) 

58.39  (56.47) 

DZP  CISD 

0.00 

93.83  (90.39) 

46.24  (44.62) 

47.59  (45.77) 

TZ2P  CISD 

0.00 

94.85  (91.54) 

46.45  (44.89) 

48.39  (46.65) 

TZ(2df,2pd)  CISD 

0.00 

95.92  (92.72) 

46.73  (45.26) 

49.19  (47.46) 

DZP  CISD+Q 

0.00 

84.60  (81.17) 

43.97  (42.35) 

40.64  (38.82) 

TZ2P  CISD+Q 

0.00 

86.24  (82.94) 

44.59  (43.03) 

41.65  (39.91) 

TZ(2df,2pd)  CISD+Q 

0.00 

87.30  (84.10) 

44.70  (43.22) 

42.60  (40.87) 

DZPCCSD 

0.00 

85.18  (81.73) 

44.80  (43.22) 

40.37  (38.51) 

TZ2P  CCSD 

0.00 

86.65  (83.22) 

45.46  (43.81) 

41.19  (39.41) 

TZ(2df,2pd)  CCSD 

0.00 

87.48  (84.15) 

45.56  (43.99) 

41.92  (40.15) 

DZPCCSD(T) 

0.00 

79.33  (75.88) 

41.95  (40.32) 

37.38  (35.56) 

TZ2P  CCSD(T) 

0.00 

79.84  (76.39) 

42.11  (40.38) 

37.73  (36.01) 

TZ(2df,2pd)  CCSD(T) 

0.00 

80.52  (77.16) 

42.17  (40.52) 

38.35  (36.64) 

a  Numbers  in  parentheses  include  zero  point  vibrational  energy  corrections.  CISD+Q 
energies  were  determined  at  the  CISD  optimized  geometries  and  corrected  for  ZPVE  using 
CISD  harmonic  vibrational  frequencies  with  the  same  basis  set. 


Chapter  IV 

Oxirene;  To  Be  or  Not  to  Be?* 


*  G.  Vacek,  J.  M.  Galbraith,  Y.  Yamaguchi,  H.  F.  Schaefer,  R.  H.  Nobes,  A.  P.  Scott  and  L.  Radom 
J.  Phys.  Chem.  98,  8660  (1994). 
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The  C2v-synmietry  structure  of  oxirene  has  been  examined  using  ab  initio 
molecular  orbital  calculations  with  large  basis  sets  and  a  variety  of  methods 
of  including  electron  correlation.  Different  qualitative  conclusions  regarding 
the  nature  of  oxirene  are  reached  depending  on  the  choice  of  basis  set  and 
method  of  electron  correlation  incorporation.  With  certain  combinations  of 
basis  set  and  theoretical  method,  the  symmetric  oxirene  structure  is  found  to 
be  a  saddle  point  on  the  potential  energy  surface,  while  for  other 
combinations  oxirene  is  a  local  minimum.  Inclusion  of  triple  excitations  in 
the  correlation  treatment  has  a  large  effect,  tending  to  make  the  curvature  of 
the  surface  corresponding  to  a  ring-opening  distortion  of  the  C2\r-symmetry 
structure  less  positive  (or  more  negative).  This  is  counterbalanced  by  a 
basis  set  effect,  with  inclusion  of/ functions  making  the  curvature  more 
positive.  At  our  highest  level  of  theory,  CCSD(T)  with  basis  sets  of 
triple-zeta  quality  and  including  multiple  d  shells  and  an/ shell  on  C  and  O 
and  multiple  p  shells  and  a  d  shell  on  H,  oxirene  is  a  genuine  minimum 
under  harmonic  vibrational  analysis,  with  a  ring-opening  frequency  of  139- 
163  cm'^ 


Introduction 

The  decomposition-rearrangement  of  diazoketones  and  diazoesters  to  form  ketenes 


and  nitrogen,  the  Wolff  rearrangement  (1),  was  first  documented  in  1902.1 
N2v 
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Although  almost  a  century  has  passed,  the  mechanism  of  the  Wolff  rearrangement  is  still  in 
doubt.  Numerous  possibilities  for  the  mechanism  have  been  proposed,  but  to  date  no 
consensus  among  chemists  has  been  reached.  Meier  and  Zeller^  have  reviewed  the 
progress  made  in  determining  the  Wolff  rearrangement  mechanism,  including  the  possible 
roles  of  both  free  and  complexed  carbenes,  1,3-dipoles,  1,3-diradicals  and  oxirenes  as 
intermediates  or  transition  structures.  Wolff  himself  believed  that  reaction  (1)  involved 
ketocarbenes  as  intermediates.  More  recently,  the  participation  of  symmetric  intermediates. 
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particularly  oxirenes,  has  been  proposed  and  supported  by  results  from  isotopic-labeling 
experiments, despite  the  fact  that  earlier  experimental  studies  failed  to  lend  support  for 
their  existence.^ 

The  history  of  the  study  of  oxirenes,  and  to  some  extent  ketocarbenes,  is  well 
documented  in  Lewars's  1983  review  ^  At  that  time,  experimental  evidence  only  supported 
the  existence  of  oxirenes  as  short-lived  intermediates  in  photochemical  and  thermal  Wolff 
rearrangements,^.^  in  the  peroxyacid  oxidations  of  acetylenes'^  and  in  ketene  photolysis 
reactions.^  Since  then,  evidence  has  been  presented  for  the  observation  of  a  substituted 
oxirene  by  laser  flash  photolysis  during  a  Wolff  rearrangement^  and  for  the  detection  by 
FT-IR  spectroscopy  of  dimethyloxirene,  isolated  in  rare-gas  matrices  during  the  photolysis 
of  3-diazo-2-butanone.io 

Several  experimental  attempts  have  been  made  to  isolate  the  parent  oxirene  molecule 
in  the  gas  phase.  Most  recently.  Hop,  Holmes  and  Terlouw^i  have  tried  to  produce 
oxirene  by  the  creation  and  subsequent  vertical  neutralization  of  oxirene  radical  cation, 
which  had  been  predicted  theoretically  to  be  stable.*^  This  attempt  was  unfortunately 
unsuccessful.  Thus,  no  firm  conclusion  as  to  whether  unsubstituted  oxirenes  are  true 
minima  on  the  potential  energy  surface  (PES)  or  whether  they  are  merely  transition 
structures  has  yet  been  reached.  Hope  for  the  isolation  of  oxirene  might  be  taken  from  the 
direct  observation  of  the  short-lived  X  vinylidene  species,  which  was  (albeit  with 
much  controversy)  predicted  to  exist  with  ab  initio  theory and  observed  by  Lineberger’s 
group  with  negative  ion  photoelectron  spectroscopy,  Similar  success  by  similar  means 
was  achieved  for  the  X  U'  state  of  monofluorovinylidene  which  has  an  isomerization 
barrier  height  of  only  2  kcal  mol-i.'^-*^  Unfortunately,  oxirene  has  no  stable  anion  and  so 
negative  ion  photoelectron  spectroscopy  will  not  be  a  viable  technique  for  observing  neutral 
oxirene. 

Numerous  theoretical  investigations,  both  semiempirical*'^-!*  and  ab  initio,^^-'^^ 
have  been  undertaken  on  the  PES  that  includes  ketene,  formylmethylene  and  oxirene. 
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MINDO/3  and  NDDO  calculations  find  oxirene  to  be  more  stable  than  formylmethylene 
by  about  20  kcal  mol-f  Ab  initio  results  at  the  Hartree-Fock  level, ‘^’22-24  on  the  other 
hand,  show  formylmethylene  to  be  more  stable  than  oxirene  by  12-17  kcal  mol-f  When 
electron  correlation  is  included  through  single-point  calculations  on  the  Hartree-Fock 
optimized  structures, 23>24  formylmethylene  is  still  favored  but  the  energy  difference  is 
reduced  to  1-8  kcal  mokf  Moreover,  such  calculations  have  suggested  that  singlet 
formylmethylene  can  collapse  without  a  barrier  to  ketene,  the  Hartree-Fock  transition 
structure  for  rearrangement  lying  lower  in  energy  on  the  correlated  surface  than 
formylmethylene  itself.23.24  We  note  that,  in  the  thermolysis  of  diazoketones,  the 
participation  of  oxirene  appears  to  be  small  at  low  temperatures  where  ketocarbenes  are 
most  important  but  increases  with  increasing  temperature.^  This  implies  that  ketocarbenes 
are  energetically  more  stable  than  oxirenes.  This  agrees  with  the  ab  initio  energy 
differences  but  casts  some  doubt  on  the  conclusion  regarding  the  stability  of 
formylmethylene  reached  on  the  basis  of  single-point  correlated  calculations.23.24 

Ketocarbenes  have  been  isolated  in  low-temperature  matrices.26  These 
ketocarbenes  were  triplet  ground  states,  however,  that  should  not  contribute  to  the  Wolff 
rearrangement  which  is  considered  to  occur  on  the  singlet  surface.22  Nevertheless,  the 
possibility  of  singlet-triplet  crossing  cannot  be  discarded  entirely.  The  interested  reader  is 
referred  to  an  MCSCF  study  by  Novoa,  McDouall  and  Robb  for  a  discussion  of  this 

possibility.28 

Those  ab  initio  studies '^"24  which  have  considered  the  barrier  to  ring  opening  for 
oxirene  have  all  concluded  that  the  barrier  is  less  than  7.3  kcal  mokf  The  most  reliable 
figure  is  perhaps  that  of  2  kcal  mob'  as  determined  by  Tanaka  and  Yoshimine  by  using  a 
reaction  coordinate  approach.23  The  existence  of  a  barrier  at  low  levels  of  theory  suggests 
that  oxirene  should  be  a  minimum,  but  a  more  rigorous  test  lies  in  vibrational  analyses  at 
higher  levels  of  theory.  Bouma  et  alP-^  determined  that  at  the  HF/3-21G  level  of  theory  the 
C2v-symmetry  form  of  oxirene  is  a  minimum  with  a  ring-opening  frequency  of  524  cm'f 
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More  recently  Vacek,  Colegrove  and  Schaefer^s  found  a  decrease  (from  524  to  445  cm-*) 
in  the  ring-deformation  (62)  vibrational  frequency  when  the  basis  set  is  improved 
moderately  (to  DZP).  The  Vacek  paper^s  also  shows  this  particular  frequency  decreasing 
further  (from  445  to  262  cm-*)  with  inclusion  of  electron  correlation.  Although  at  their 
highest  level  of  theory  (CCSD/DZP)  they  still  conclude  that  oxirene  is  a  minimum,  Vacek 
et  al.'^  were  not  confident  that  the  effects  of  still  larger  basis  sets  and  improved  treatments 
of  electron  correlation  would  not  cause  this  particular  frequency  to  become  imaginary. 

Clearly,  there  are  many  questions  that  remain  unanswered  concerning  oxirene  and 
related  species  that  would  benefit  from  a  high  level  ab  initio  determination  of  the  entire 
oxirene  — »  ketene  potential  energy  hypersurface,  and  from  calculations  on  related  larger 
homologues.  These  include  the  issue  of  substitution  effects  on  the  stability  of  oxirenes,^^ 
an  understanding  of  the  kinetics  of  reactions  on  the  C2H2O  PES^O-^*  and  further  probing 
of  the  mechanism  for  the  Wolff  rearrangement.32  fo  the  present  paper,  we  address  the 
particular  aspect  of  whether  or  not  high  level  ab  initio  calculations  indicate  oxirene  to  be  a 
true  minimum  on  the  surface. 

Theoretical  Methods 

The  C2v~symmetry  stmcture  of  oxirene  was  studied  with  a  variety  of  basis  sets  and 
methods  of  electron  correlation.  Calculations  were  carried  out  using  the  Gaussian 
92/DFT33  and  ACES  1134  programs  (at  the  Australian  National  University,  ANU)  and  the 
PSI  program  (at  the  Center  for  Computational  Quantum  Chemistry,  CCQC).35 

Calculations  at  the  ANU  employed  the  standard  6-31G(d),  6-31G(J,p)  and 
6-311G(d,p)  basis  sets36and  the  Dunning  correlation-consistent  basis,  cc-pVTZ.37 
Additional  polarization  functions  and  diffuse  functions  were  added  to  the  6-31G(d,p)  and 
6-311G(<f,p)  basis  sets  to  form  the  6-Z\G{df,p),  6-311G(4/’,p)  and  6-311+G(2(i/,p) 
sets.38.39  Basis  sets  formed  from  6-3 IG  employed  sets  of  six  Cartesian  d-like  and  ten 
Cartesian /-like  polarization  functions  while  those  based  on  6-31 IG  and  cc-pVTZ 


68 


employed  five  pure  d  and  seven  pure/polarization  functions,  i.e.  spherical  harmonics. 

The  DZP  basis  used  in  the  calculations  performed  at  the  CCQC  was  derived  from 
the  standard  DZ  basis  by  augmenting  with  a  set  of  d-Viks,  functions  on  the  heavy  atoms  and 
a  set  of  p  functions  on  hydrogen.38.40  The  DZP++  set  was  constructed  by  adding  diffuse  ^ 
and/?  functions  to  the  DZP  basis.4i-42  The  TZ2P  basis  was  constructed  from  the  standard 
TZ  basis  by  the  addition  of  two  sets  of  polarization  functions  on  each  atom.38.43  This 
basis  was  further  enhanced  by  the  addition  of /-like  polarization  functions  on  the  heavy 
atoms  and  a  set  of  rf-like  functions  on  hydrogen, 38.44  leading  to  the  basis  referred  to  as 
TZ2?(f,d)  in  the  text.  Yet  another  basis,  TZ2P++,  was  formed  by  the  addition  of  diffuse 
functions  to  all  atoms  in  the  TZ2P  set.4i.45  The  QZ3?(f,d)  basis  was  constructed  by  the 
addition  of  higher  angular  momentum  functions  to  a  standard  QZ  basis.38.46  All  the  CCQC 
basis  sets  employed  sets  of  six  Cartesian  rf-like  and  ten  Cartesian  /-like  polarization 
functions. 

The  molecular  structure  of  oxirene  was  fully  optimized  using  analytic  gradient 
techniques  for  restricted  Hartree-Fock  (RHF),47  two-configuration  SCF  (TCSCF),48 
M0ller-Plesset  perturbation  theory  (MPn)49-5i,  configuration  interaction  (CI),32.53 
coupled-cluster  (CC)34.55  and  quadratic  configuration  interaction  (QCI)36.57 
wavefunctions.  Density  functional  calculations  were  also  performed  using  the  local  SVWN 
functional,  38  the  gradient-corrected  BLYP  functional  38.59  and  the  hybrid  BeckeSLYP 
functional.^!!  In  all  cases,  the  default  geometry  optimization  criteria  were  employed.  For 
the  iterative  methods,  all  single  and  double  excitations  (SD)  from  the  SCF  reference 
configuration  were  included  (CISD,  CCSD  and  QCISD).  For  the  coupled-cluster  and 
quadratic  configuration  methods,  the  effects  of  connected  triple  excitations  were  included 
perturbatively  [CCSD(T),  QCISD(T)].35,57  in  niost  of  the  calculations  performed  at  the 
ANU,  all  orbitals  were  correlated  (full,  abbreviated  fu),  whereas  in  the  CCQC  calculations 
only  the  valence  electrons  were  explicitly  correlated,  i.e.  the  three  lowest-lying  MOs  (C  and 
O  Is-like)  were  held  doubly  occupied  (frozen  core,  abbreviated  fc). 
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For  the  TCSCP  wavefunctions,  the  second  configuration  was  a  double  excitation 
from  the  highest  occupied  MO  {2b\)  to  the  lowest  unoccupied  MO  (102)-  This  is  a  ;r-> 

excitation. 

The  RHF,  density  functional  and  MP2  quadratic  force  constants  were  evaluated  via 
analytic  second  derivative  procedures. The  method  of  finite  central  differences  of  the 
analytically  derived  gradients  was  utilized  to  evaluate  the  quadratic  force  constants  for  the 
remaining  wavefunctions  used  in  this  study. 

Results  and  Discussion 

Structural  Predictions 

Total  energies  and  structural  predictions  are  shown  in  Table  I.  Perhaps  the  most 
interesting  feature  is  the  relative  insensitivity  of  the  geometrical  parameters  to  basis  set 
improvements  for  a  given  theoretical  method.  For  example,  with  the  RHF  method,  the 
C-O  bond  length  falls  within  a  range  of  0.01  A  across  the  entire  set  of  bases.  This 
contrasts  with  the  large  change  seen  in  going  from  the  smaller  basis  set  4-3 IG  results  (rco 
=  1.552  A)  of  Tanaka  and  Yoshimine  to  their  DZ  +  P  results  (rco  =  1-465  A),  where  there 
is  a  change  of  0.09  A.23  Clearly  the  first  set  of  polarization  functions  is  very  important. 
On  the  other  hand,  increasing  the  size  of  the  sp  part  of  the  basis  (from  a  double  to  a  triple 
split  of  the  valence  functions  or  by  adding  a  set  of  diffuse  functions)  has  very  little  effect  on 
the  calculated  geometry.  However,  the  addition  of /-like  polarization  functions  does  have  a 
significant  effect,  especially  on  the  C-O  bond  length  (decreasing  this  length  by  0.010- 
0.015  A  for  the  conventional  correlated  methods).  It  is  somewhat  reassuring  to  see  the 
geometrical  results  converge  so  quickly  with  basis  set.  Unfortunately,  as  discussed  below, 
the  harmonic  vibrational  frequencies  are  not  so  convergent  with  respect  to  increase  in  basis 
set  size. 

The  changes  in  geometrical  parameters  with  increasing  sophistication  of  theoretical 
method  are  somewhat  greater  than  with  basis  set.  However,  the  trends  observed  (namely 
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an  increase  in  bond  lengths  as  one  improves  the  electron  correlation  scheme  ^7,68)  are 
typical  for  C-0,  C-C  and  C-H  bond  distances.  For  the  correlated  methods,  results  with 
all  orbitals  included  (fu)  are  very  close  to  those  using  the  frozen-core  approach  (fc)  in  the 
few  cases  for  which  direct  comparisons  are  possible.  The  inclusion  of  triple  excitations  in 
the  correlation  treatment  [MP4SDQ  — >  MP4,  QCISD  QCISD(T),  CCSD  >  CCSD(T)] 

has  a  relatively  large  effect,  tending  to  increase  the  lengths  of  the  C-C  and  C-O  bonds. 

0  . 

Note  also  the  large  fluctuations  in  the  Mpller-Plesset  series,  with  a  change  of  0.021  A  in 
the  C-0  length  in  going  from  MP3  to  MP4,  and  the  wide  range  of  values  for  the  C-O 

length  from  the  density  functional  approaches. 

Despite  the  relatively  large  changes  in  the  C-O  lengths,  it  appears  that  the 
predictions  of  geometrical  structure  are  converging  towards  a  single  set  of  values  with 
electron  correlation.  The  largest  part  of  the  change  in  geometry  occurs  between  the  SCF 
and  either  the  MP2  or  the  CISD  methods  of  electron  correlation.  More  sophisticated 
inclusion  of  electron  correlation  has  a  much  smaller  effect. 

Frequency  Predictions 

Predictions  of  the  harmonic  vibrational  frequencies  appear  in  Table  H.  Of  particular 
interest  is  the  ring-deformation  normal  vibrational  mode.  It  can  be  seen  that  the  nature  of 
the  oxirene  stationary  point  changes  depending  on  choice  of  basis  set  and  theoretical 
procedure'.  If  one  were  to  use  only  SCF  calculations,  the  clear  conclusion  would  be  that 
oxirene  is  a  minimum  on  the  surface.  Conversely,  all  of  the  density  functional  based 
approaches  indicate  that  oxirene  is  a  saddle  point.  For  the  conventional  correlated 
procedures,  one  sees  a  mixture  of  real  and  imaginary  values  for  the  ring  deformation 
vibrational  mode. 

We  note  initially  that  the  effect  of  freezing  the  core  in  the  correlated  calculations  is 
very  small.  For  example,  the  mean  absolute  difference  between  the  QCISD-fc/6-31G(if) 
and  QCISD-fu/6-31G(£0  frequencies  is  just  3  cm*!,  with  a  largest  difference  of  13  cm-1. 
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The  drop  in  the  ring-opening  harmonic  vibrational  frequency  from  RHF/DZP  (445 
cm-i)  to  CISD/DZP  (338  cm-*)  to  CCSD/DZP  (262  cm-*)  has  already  been  noted  by  Vacek 
et  al.^  Inclusion  of  triple  excitations  in  the  correlation  treatment  has  a  large  effect,  further 
flattening  the  surface  at  the  C2v-synimetry  stationary  point.  For  example,  for  the  DZP 
basis  set,  there  is  another  drop  (to  119  cm-*)  when  going  to  the  CCSD(T)  method  of 
electron  correlation.  Furthermore,  when  making  rather  typical  basis  set  improvements 
[such  as  from  DZP  to  TZ2P  or  6-3lG(d,p)  to  6-31  lG(d,p)],  we  again  see  a  lowering  of 
this  vibrational  frequency.  This  effect  is  somewhat  additive  and  as  a  result  it  can  be  seen 
that  at  the  CCSD(T)  level  with  large  spd  basis  sets  the  surface  actually  has  negative 
curvature,  i.e.  oxirene  is  a  saddle  point  at  these  levels  of  theory! 

However  all  other  basis  set  improvements  seem  to  stabilize  oxirene  with  respect  to 
ring  deformation.  The  addition  of  diffuse  functions  to  the  DZP  and  TZ2P  bases  slightly 
decreases  the  C-O  bond  length  and  stabilizes  the  mode  in  question  by  12  and  6  cm-*, 
respectively,  on  the  RHF  surfaces  and  by  56  and  72  cm-*  on  the  CCSD(T)  surfaces.  Note 
also  that,  although  the  initial  valence  and  polarization  improvement  from  DZP  to  TZ2P  was 
destabilizing,  an  additional  improvement  from  TZ2P(f,d)  to  QZ3P(f,d)  is  stabilizing,  at 
least  on  the  RHF  surface. 

The  stabilization  of  oxirene  by  higher  angular  momentum  functions  is  even  greater 
than  that  by  diffuse  functions.  For  example,  with  the  CCSD(T)  method,  the  addition  of  an 
/shell  to  the  6-311G(d,p)  basis  or/ (on  C  and  O)  and  d  (on  H)  shells  to  the  TZ2P  basis 
has  a  large  effect,  restoring  positive  curvature  to  the  potential  surface.  With  our  best 
theoretical  levels,  CCSD(T)-fc/TZ2P(/;iO  and  CCSD(T)-fu/cc-pVTZ,  oxirene  is  a  local 
minimum  on  the  potential  surface  with  ring-opening  frequencies  of  163  and  139  cm-*, 
respectively. 

The  poor  performance  of  Mpller-Plesset  perturbation  theory  for  this  problem 
should  be  noted.  One  can  see  rather  large  fluctuations  in  the  b2  vibrational  frequency  along 
the  series  RHF  ^  MP2  -»  MP3  ^  MP4  (459  135  ^  400  -4  207  i  cm-*)  for  the 
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6-3 1  \G{df,p)  basis.  The  importance  of  a  balanced  treatment  of  triple  substitutions  is  also 
nicely  illustrated  by  results  with  the  6-2\\G{df,p)  basis  set.  MP4  without  the  triples 
(MP4SDQ)  has  a  ring-opening  frequency  of  227  cm'^  This  changes  dramatically  to  207  i 
cm'^  when  triples  are  included  (MP4).  MP4  is,  however,  known^^70  jq  overestimate  the 
importance  of  the  triple  substitutions.  A  more  complete  treatment  of  the  triples  is  provided 
by  QCISD(T)  for  which  the  ring-opening  frequency  becomes  real  (82  cm-l).  CCSD(T) 
results  in  a  further  improvement  in  the  treatment  of  the  triple  effects, and  at  the 
CCSD(T)/6-31  lG(cl/,p)  level  the  bi  frequency  increases  to  137  cm'*.  The  effect  of  the 
inclusion  of  quadruple  excitations^i  on  the  geometry  and  the  ring-opening  frequency  of 
oxirene  would  certainly  be  of  interest! 

A  brief  aside  should  be  made  concerning  the  possible  importance  of  a 
multireference  wavefunction.  For  all  three  CISD  wavefimctions  constructed  in  this  study, 
the  coefficient  of  the  reference  configuration  is  0.94.  The  coefficients  for  the  second  most 
important  configuration  (a  double  excitation  from  the  2b\  ;it;c-bonding  MO  to  the  1^2 
;r*cc-antibonding  MO)  are  -0.06,  -0.05  and  -0.04  for  increasing-sized  basis  sets.  Such 
contributions  from  the  second  configuration  imply  that  single-reference  CISD  may  not  be 
entirely  suitable.  However,  analysis  of  the  coupled-cluster  wavefunctions  reveal  T\ 
diagnostic  values  of  0.014  for  all  cases.  This  is  less  than  the  suggested  threshold  of 
0.02, "^2  implying  that  most  of  the  dynamical  correlation  energy  is  being  recovered  by  the 
coupled-cluster  wavefunctions.  Thus  while  CISD  results  are  suggested  to  be  insufficient, 
coupled-cluster  wavefunctions  should  be  adequate  to  resolve  the  problem  without  resorting 
to  multireference  techniques.  This  last  conclusion,  of  course,  only  pertains  to  the  oxirene 
stationary  point  and  does  not  necessarily  apply  to  the  entire  isomerization  hypersurface. 

Conclusion 

The  present  study  has  examined  the  C2v-synimetry  structure  of  oxirene  with  a 
variety  of  basis  sets  and  theoretical  procedures.  The  results  at  the  various  levels  show 
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qualitative  differences,  specifically  with  regard  to  the  nature  of  the  stationary  point. 
Density  functional  methods,  fourth-order  Mpller-Plesset  perturbation  theory  and  CCSD(T) 
methods  with  large  spd  basis  sets  all  show  the  symmetric  structure  to  be  a  transition 
structure  along  the  ring-opening  normal  coordinate,  while  other  methods  show  it  to  be  a 
true  minimum.  At  our  highest  level  of  theory,  CCSD(T)  with  basis  sets  of  triple-zeta 
quality  and  including  multiple  d  shells  and  an  f  shell  on  C  and  O  and  multiple  p  shells  and  a 
d  shell  on  H,  oxirene  is  a  true  minimum  with  a  ring-opening  frequency  of  139-163  cm-l. 

The  large  basis  set  CCSD(T)  calculations  may  prove  prohibitively  costly  for 
studying  the  low-symmetry  structures  in  the  oxirene  ketene  isomerization.  A 
compromise  in  terms  of  computational  efficiency  and  accuracy  may  be  provided  by  the 
CCSD(T)/6-311G(rf/,/7)  level  of  theory,  which  predicts  a  geometry  and  a  ring-opening 
frequency  satisfactorily  close  to  those  given  by  the  CCSD(T)/cc-pVTZ  and 
CCSD(T)/TZ2P(/;cO  methods. 
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TABLE  I.  Calculated  total  energies  and  molecular  structures  for  the  C2v-symmetry 
structure  of  oxirene.^ 


o 


Level  of  theory*’’^ 

Src  d  Tot.  E.e 

r(C  -  0) 

r(C  =  C) 

r(C  -  H) 

<CCH 

RHF/6-31G(J) 

ANU  -1.58390 

1.467 

1.244 

1.062 

162.8 

RHF/DZP 

CCQC  -1.61718 

1.466 

1.251 

1.066 

162.6 

RHF/DZP++ 

CCQC  -1.62189 

1.464 

1.251 

1.066 

162.7 

RHF/6-311G(rf,p) 

ANU  -1.62365 

1.464 

1.243 

1.061 

162.8 

RHF/6-311G(4/ip) 

ANU  -1.63164 

1.459 

1.241 

1.061 

162.6 

RHF/TZ2P 

CCQC  -1.64283 

1.465 

1.240 

1.059 

162.2 

RHF/TZ2P++ 

CCQC  -1.64395 

1.464 

1.240 

1.059 

162.3 

RHF/cc-pVTZf 

ANU  -1.64536 

1.460 

1.240 

1.059 

162.5 

RHF/TZ2P(/:if) 

CCQC  -1.65069 

1.459 

1.240 

1.059 

162.5 

RHF/QZ3P(/:cO 

CCQC  -1.65686 

1.457 

1.240 

1.059 

162.5 

TCSCF/DZP 

CCQC  -1.64882 

1.461 

1.274 

1.065 

162.8 

TCSCF/TZ2P 

CCQC  -1.67329 

1.460 

1.261 

1.058 

162.5 

SVWN/6-311G(4f,p) 

ANU  -1.69702 

1.481 

1.265 

1.081 

161.8 

BLYP/6-311G(rf,p) 

ANU  -2.47827 

1.530 

1.273 

1.076 

161.7 

BLYP/6-311G(d/p) 

ANU  -2.48410 

1.525 

1.271 

1.077 

161.4 
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BLY?/6-3n+G{2df.p) 

Becke3LYP/6-311G(4^/7) 

MP2-fc/6-31G(^ 

MP2-fc/6-311G(i/,p) 

MP2-fc/6-311G(4/;p) 

MP2-fu/6-311G(4/;;7) 

MP2-fc/6-311+G(2ti//?) 

MP2-fu/cc-pVTZf 

MP3-fu/6-311G(d/p) 

MP4SDQ-fu/6-3 1  lGidf,p) 
MP4-fu/6-311G(4/:p) 

CISD-fc/DZP 

CISD-fc/TZ2P 

CISD-fc/TZ2P(/;c/) 

QCISD-fc/6-31G(<i) 

QCISD-fu/6-31G(ti) 

QCISD-fu/6-311G(J/p) 

QCISD(T)-fc/6-31G(iO 
QCISD(T)-fui/6-31G(rf) 
QCISD(T)-fu/6-3 1 1  Gidf,p) 


ANU 

-2.49660 

1.523 

ANU 

-2.51811 

1.499 

ANU 

-2.01868 

1.512 

ANU 

-2.09517 

1.504 

ANU 

-2.14680 

1.492 

ANU 

-2.20418 

1.489 

ANU 

-2.18535 

1.500 

ANU 

-2.23792 

1.497 

ANU 

-2.21512 

1.475 

ANU 

-2.22060 

1.482 

ANU 

-2.24683 

1.496 

CCQC 

-2.02266 

1.491 

CCQC 

-2.09860 

1.488 

CCQC 

-2.14890 

1.475 

ANU 

-2.03619 

1.506 

ANU 

-2.04774 

1.504 

ANU 

-2.22152 

1.483 

ANU 

-2.05324 

1.514 

ANU 

-2.06494 

1.512 

ANU 

-2.24511 

1.490 

1.270 

1.077 

161.6 

1.261 

1.071 

161.7 

1.277 

1.074 

161.9 

1.277 

1.071 

162.2 

1.274 

1.073 

161.9 

1.272 

1.072 

161.8 

1.271 

1.071 

162.0 

1.268 

1.064 

161.8 

1.266 

1.071 

161.8 

1.268 

1.072 

161.9 

1.277 

1:074 

161.9 

1.272 

1.073 

162.0 

1.255 

1.062 

161.5 

1.253 

1.061 

162.1 

1.275 

1.075 

162.1 

1.273 

1.075 

162.2 

1.269 

1.072 

162.0 

1.281 

1.077 

162.0 

1.280 

1.077 

161.9 

1.276 

1.074 

161.9 
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CCSD-fc/DZP 

CCSD-fu/6-311G(4/;p) 

CCSD-fc/TZ2P 

CCSD(T)-fu/6-31G(40 

CCSD(T)-fu/6-31GWp) 

CCSD(T)-fu/6-3 1  G{df,p) 

CCSD(T)-fc/DZP 

CCSD(T)-fc/DZP++ 

CCSD(T)-fu/6-3 1  \G{d,p) 

CCSD(T)-fu/6-3 1  \G{df,p) 

CCSD(T)-fc/TZ2P 

CCSD(T)-fc/TZ2P-H- 

CCSD(T)-fu/cc-pVTZf 

CCSD(T)-fc/TZ2P(/,if) 


CCQC  -2.07317 

1.506 

ANU  -2.22005 

1.481 

CCQC  -2.15573 

1.505 

ANU  -2.06428 

1.512 

ANU  -2.08064 

1.510 

ANU  -2.16296 

1.500 

CCQC  -2.09085 

1.515 

CCQC  -2.10026 

1.513 

ANU  -2.18937 

1.504 

ANU  -2.24439 

1.490 

CCQC  -2.18130 

1.518 

CC(3C  -2.18349 

1.517 

ANU  -2.27586 

1.497 

CCQC  -2.23948 

1.503 

1.285 

1.078 

161.8 

1.269 

1.072 

161.9 

1.267 

1.067 

161.3 

1.279 

1.076 

161.9 

1.280 

1.070 

161.6 

1.269 

1.071 

161.8 

1.292 

1.080 

161.7 

1.291 

1.079 

161.9 

1.280 

1.073 

161.9 

1.276 

1.074 

161.8 

1.275 

1.069 

161.2 

1.275 

1.069 

161.2 

1.270 

1.065 

161.9 

1.273 

1.069 

161.8 

a  Energies  in  hartrees,  bond  lengths  in  angstroms,  bond  angles  in  degrees. 

^  The  notation  “fu”  (full)  indicates  that  all  orbitals  were  included  in  the  correlation 
treatment  while  “fc”  (frozen-core)  indicates  that  the  core  orbitals  were  held  frozen, 
c  Calculations  based  on  the  6-3 1 G  or  CCQC  basis  sets  used  6d  and/or  1 0/ while 
calculations  based  on  6-31 IG  or  cc-pVTZ  used  5d  and/or  If. 
d  ANU;  Australian  National  University,  Canberra;  CCQC:  Center  for  Computational 
Quantum  Chemistry,  Athens. 

®  Total  energy  -150  hartrees. 
f  {As3p2d\f\ 


TABLE  II.  Calculated  harmonic  vibrational  frequencies  for  the  C2v-symmetry  structure  of  oxirene.' 
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Chapter  V:  Concluding  Remarks 


When  I  started  doing  theoretical  chemistry  in  the  summer  of  1990, 1  thought  that  the 
computer  would  do  all  the  work  for  me.  Here  I  was  only  partly  right.  The  computer  can 
solve  the  equations  that  I  tell  it  to  (within  certain  limits,  of  course),  as  far  as  handling  the 
integrals  and  summations  and  what  not.  However,  it  is  not  as  simple  as  all  that. 

Firstly,  a  project  needs  to  be  chosen  for  which  theoretical  methods  can  reasonably 
be  expected  to  add  sufficient  insight  to  illuminate  the  chemical  problem.  While  chemical 
problems  are  numerous,  not  all  fit  the  above  standards.  For  instance,  many  chemical 
systems  are  too  large  to  be  accurately  studied  with  current  computer  resources.  Fortunately 
for  those  of  us  in  this  field,  computers  are  getting  "larger"  at  a  very  good  rate.  Projects 
which  were  difficult  to  study  on  the  campus  mainframe  when  I  started  five  years  ago  can 
now  be  handled  on  a  personal  computer.  Even  if  the  chemical  system  is  small  enough, 
standard  methods  will  have  somewhat  more  difficulty  with  excited  states  and  weakly  bound 
complexes  and  certain  other  systems. 

Provided  with  a  suitable  chemical  system,  one  faces  the  second  challenge  of 
choosing  the  best  level  of  theory  with  which  to  study  the  system.  Generally  the  trade-off  is 
between  relative  accuracy  of  results  and  consumption  of  computer  resources.  That  is,  low 
accuracy  with  little  consumption  or  higher  accuracy  at  greater  expense.  To  make  matters 
worse,  there  is  no  a  priori  way  to  know  how  accurate  a  particular  method  will  be  for  a 
particular  system.  Experience  and/or  systematic  studies  on  related  systems  will  help  by 
providing  guidelines  to  this  choice,  but  often  one  must  just  keep  beating  the  problem  with 
successively  bigger  theoretical  hammers  until  the  desired  degree  of  accuracy  is  achieved. 
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The  final  consideration  in  computational  theoretical  chemistry  which  prevents  the 
computer  from  completing  the  research  for  us  is  the  nature  of  the  input  and  output  data.  As 
with  any  computation,  garbage  in  means  garbage  out.  A  decent  chemical  understanding  of 
the  system  of  interest  is  therefor  a  must  to  give  the  correct  input.  Furthermore,  chemical 
intuition  is  necessary  to  sort  the  output  data  in  some  meaningful  manner. 

So,  as  I  stated,  theoretical  computational  chemistry  turns  out  to  not  be  quite  so 
simple  as  I  had  believed  prior  to  starting  in  the  field.  However,  the  fact  that  it  is  no  so  easy 
is  the  only  thing  which  keeps  it  from  getting  boring  (and  me  from  getting  replaced  by  said 
computer).  Besides,  I'll  take  this  over  longhand  calculations  on  the  chalkboard  any  day. 

I  will  end  with  some  closing  remarks  about  the  individual  projects  that  comprise 
this  work. 

Within  my  earshot,  acetylene  (really  all  C2H2  systems)  was  referred  to  as  "the  little 
system  that  could."  These  systems  have  been  studied  by  chemists  since  the  dawn  of 
chemistry,  and  more  recently  by  theoretical  chemists  since  the  dawn  of  theoretical 
chemistry.  Yet  there  still  remain  many  interesting  questions  to  be  answered  about  these 
systems,  not  to  mention  the  many  related  A2H2  systems  for  which  C2H2  is  a  model.  Here 
we  have  demonstrated  that  for  these  systems  in  order  to  obtain  relative  energetic  predictions 
to  within  an  accuracy  of  2  kcal/mole,  it  is  necessary  to  use  basis  sets  with  multiple 
polarization  functions  as  well  as  higher-angular-momentum  polarization  functions  and 
include  electron  correlation  from  higher  than  singles  and  doubles  excitations.  The 
isomerization  reaction  from  b  ^Bu  trans-hexA  to  a  ^^2  cw-bent  acetylene  is  predicted  to  be 
exothermic  by  8.3  kcal/mole  with  a  barrier  height  of  only  4.8  kcal/mole.  This  may  help 
explain  why  the  trans  state  has  yet  to  be  experimentally  characterized.  Hopefully  our 
definitive  theoretical  predictions  will  provide  a  stimulus  for  the  full  experimental 
characterization  of  the  lowest  triplet  states  of  C2H2. 

All  the  research  in  this  theoretical  work  has  been  completed  with  the  thought  in 
mind  that  it  would  influence  and  aid  not  only  other  theoretical  chemists  but  also 
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experimental  chemists.  Thus  it  is  all  the  more  gratifying  that  the  study  of  AlOH-HAlO 
helped  experimentalists  whom  I  know  personally.  Our  predictions  compared  sufficiently 
well  to  their  experimental  results  so  as  to  confirm  the  assignments  which  they  determined. 
Not  all  projects  work  out  this  fmitfully. 

The  oxirene  study  was  just  plain  nasty.  Results  at  various  levels  of  theory  show 
differences  not  just  of  quantitative  accuracy  but  of  the  qualitative  nature  of  the  stationary 
point.  At  the  highest  levels  of  theory  oxirene  is  shown  to  be  a  true  minimum. 
Experimental  confirmation  of  this  hard  won  fact  would  be  wonderful.  Matrix-isolation 
infrared  (IR)  spectroscopy  is  an  obvious  choice  for  the  laboratory  identification  of  oxirene. 
However,  observation  by  this  technique  may  be  problematic.  For  example,  the  vinylidene 
species  with  a  barrier  to  isomerization  of  roughly  2  kcal/mole  is  well  established  as  a 
molecular  entity  by  negative-ion  photodetatchment  spectroscopy  but  has  yet  to  be  observed 
via  matrix-isolation  IR  spectroscopy.  Oxirene  has  no  stable  negative-ion,  as  well,  so 
photodetatchment  spectroscopy  is  not  an  option. 

The  high  levels  of  theory  used  in  this  work  to  definitively  show  oxirene  to  be  a 
minimum  are  prohibitively  costly  for  the  study  of  the  other  (lower-symmetry)  structures  on 
the  oxirene-ketene  isomerization  reaction  surface  or  the  Wolff-rearrangement  reaction 
surface.  We  have  suggested  compromise  levels  of  theory  which  are  both  accurate  enough 
and  efficient  enough  for  that  purpose.  Such  levels  of  theory  have  been  taken  up  by  others 
and  lead  to  fruitful  studies  of  the  entire  PE  hypersurface.'  Those  results  clearly  support  the 
participation  of  the  oxirene  minimum  in  the  Wolff  rearrangement.  Entry  into  the 
rearrangement  via  formylmethylene  will  still  result  in  carbon  scrambling  via  oxirene  since 
the  energetic  barrier  to  scrambling  (oxygen-ring  formation  and  subsequent  opening)  is 
much  smallex  that  the  5.5  kcal/mole  barrier  to  (hydrogen  transfer)  isomerization  to  ketene. 


'  A.  P.  Scott,  R.  H.  Nobes,  H.  F.  Schaefer  and  L.  Radom,  J.  Am.  Chem.  Soc.  116,  10165  (1994). 


